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ABSTRACT 

We present evidence for spatially extended low surface brightness emission around Lyman break galaxies 
(LBGs) in the V band image of the Hubble Ultra Deep Field, corresponding to the z ~ 3 rest-frame FUV light, 
which is a sensitive measure of star formation rates (SFRs). We find that the covering fraction of molecular gas 
at z ~ 3 is not adequate to explain the emission in the outskirts of LBGs, while the covering fraction of neutral 
atomic-dominated hydrogen gas at high redshift is sufficient. We develop a theoretical framework to connect 
this emission around LBGs to the expected emission from neutral H I gas i.e., damped Lya systems (DLAs), 
using the Kennicutt-Schmidt (KS) relation. Working under the hypothesis that the observed FUV emission 
in the outskirts of LBGs is from in situ star formation in atomic-dominated hydrogen gas, the results suggest 
that the SFR efficiency in such gas at z ~ 3 is between factors of 10 and 50 lower than predictions based on 
the local KS relation. The total SFR density in atomic-dominated gas at z ~ 3 is constrained to be ~ 10% of 
that observed from the inner regions of LBGs. In addition, the metals produced by in situ star formation in the 
outskirts of LBGs yield metallicities comparable to those of DLAs, which is a possible solution to the "Missing 
Metals" problem for DLAs. Finally, the atomic-dominated gas in the outskirts of galaxies at both high and low 
redshifts has similar reduced SFR efficiencies and is consistent with the same power law. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — galaxies: 
photometry — general: galaxies — quasars: absorption lines 



1. INTRODUCTION 

Understanding how stars form from gas is vital to our com- 
prehension of galaxy formation and evolution. Although the 
physics involved in this process is not fully understood, we 
do know something about the principal sites where star for- 
mation occurs and where the gas resides. Most of the known 
star formation at high redshift occurs in Lyman break galax- 
ies (LBGs), a population of star-forming galaxies selected for 
their opacity at the Lyman limit and the presence of upper 
main sequence stars that emit FUV radiation. They also have 
very high star formation rates (SFRs) of ~8 M Q yr" 1 after 
correcting for extinction (Shapley et al. 2003t). 

While LBGs have a wide range of morphologies, the aver- 
age half-light radius for z ~ 3 LBGs is about ~ 2— 3 kpc in the 
optica l at ~ 25 mag (e.g., iGiavalisco et al.|[T996t lLaw et all 
2007). However, studies of the Hubble Ultra Deep Field 
(UDF) have shown that fainter LBGs have smaller half-light 
radii, around ~ 1 kpc for LBGs with brightn esses similar to 
those used in this study (V ~ 26 - 27 mag; iBouwens et al.1 
2004). While we have no empirical knowledge about star for- 
mation in the outer regions of LBGs, sim ulations suggest that 
stars may be forming further out (e.g.. iGnedin & Kravtsovl 
l2010h . It still remains an unanswered question whether star 
formation occurs in the outer disks of high redshift galaxies. 
Local galaxies at z ~ are forming stars in their outer disks, 
as obser ved in the ultraviolet (Thilker et al. 2005; Bigi el et aT] 
1201 Obi ah . At low redshift, thi s star formation occurs in 
atomic-dominated hydrogen gas (Fuma galli & Gavazzill2008l; 
Bigiel et al. 2010b a), where the majority of the hydrogen gas 
is atomic but molecules are present. At high redshift such gas 
resides in damped Lya systems (DLAs). 



DLAs are a population of H I layers selected for their neu- 
tral hydrogen column densities of A/hi > 2 x 10 20 cm" 2 , which 
dominate the neutral-gas content of the universe in the redshift 
interval < z < 5. In fact, DLAs at z ~ 3 contain enough gas 
to account for 25%-50% of the mass conte nt of visible matter 
in modern galaxies (see W olfe et al.l 120051 for a review) and 
are neutral-gas reservoirs for star formation. 

The locally established Ken nicutt-Schmidt (KS) relation 
dKennicuttl 1 998at ISchmidtl 1 9591) relates the SFR per unit area 
and the total gas surface density (atomic and molecular), Ssfr 
oc Sg^. While it is reasonable to use this relationship at low 
redshift in normal star-forming galaxies, many cosmological 
simulations use it at all redshifts with out distinguishing be- 
tween atomic and molecular gas (e.g.,[N agamine et alj|2004 , 
20071 120101 IRazoumov et al.ll2006t iBrooks et alj|2007l 2009 ; 
Pontzen et alJl2008t IRazoumov et alj|2008t iRazoumovl 2009 ; 
Tescari et al.l l2009t IDekel et al.l I2009albi keres et al l 12009 : 
iBarnes & Haehneltll^oi oPl. Yet at z ~ 3, excluding regions 
immediately surrounding high surface brightness LBGs, the 
SFR per unit comoving volume, p'*, of DLAs was found to 
be less than 5% of what is expected from the KS relation 
(IWolfe&Chenl 120061) . This means that a lower level of in 
situ star formation occurs in atomic-dominated hydrogen gas 
at z ~ 3 than in modern galaxiefl 

These results have multiple implications affecting such gas 
at high redshift. First, the lower SFR efficiencies in DLAs 
are inconsistent with the 158/im cooling rates of DLAs with 
purely in situ star f ormation. Specifical ly, I Wolfe etafl (120081) 
adopt the model of IWolfe et al.l (l2003bl) . in which star forma- 

3 We note that there are also a large number of papers that do not as- 
sume the KS relation in their simulations and models (e.g.,lKravtsov|[200: 
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Krumholz et al. 2008, 2009a b; Tassis et al. 2008; Robertson & Kravtsov 
l2008tlGnedin & Kravtsov<20*10Tl201 lUFeldmann et alJ201 II) . 

4 We note that the SFR efficiency discussed in this paper relates to the nor- 
malization of the KS relation, and is not the same as the star formation effi- 
ciency (SFE), which is the inverse of the gas depletion time (e.g., L erov et al] 
2008). 
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tion generates FUV radiation that heats the gas by the grain 
photoelectric mechanism. Assuming thermal balance, they 
equate the heating rates to the [C II] 158 /im cooling rates of 
DLAs i nferred from the me asured C IPA1335.7 absorption of 
DLAs dWorfeetal.ll2003bl). The DLA c ooling rates exhibit 
a bimodal distribution (IWolfe et al.ll2008h . and the population 
of DLAs with high cooling rates have inferred heating rates 
significantly higher than that implied by the upper limits of 
FUV emission of spatially extended sources (IWolfe & C hen 
120061: IWolfe et all 120081) . Therefore, another source of heat 
input is required, such as compact star-forming regions em- 
bedded in the neutral gas; e.g., LBGs. Second, since is 
directly proportional to the metal production rate, the limits 
on p„ shift the problem of metal overproduction in DLAs by 
a factor of 10 dPettinill 1999L l200l 120061: IWolfe et alJl2003bl 
; known as the "Missing Metals" probl em for DLAs), to on e 
of underproduction by a factor of three (Wol fe & Chenl2006h . 
Lastly , the multi-component ve locity structure of the DLA gas 
(e.g JProchaska &~W olfe 1997) suggests that energy input by 
supernova explosions is required to replenish the turbulent ki- 
netic energy lost through cloud collisions; i.e., some in situ 
star formation should be present in DLAs. 

One possible way to reconcile the lack of detected in situ 
star formation in DLAs at high redshift and the properties 
of DLAs that require heating of the gas, is that compact 
LBG cores embedded in spatially extended DLA gas may 
cause both the heat input, chemical enrichment, and turbu- 
lent kinetic energy observed in DLAs. There are several in- 
dependent lines of evidence linking DLAs and LBGs; e.g. 
(1) there is a significant cros s correlation between LBGs 
and DLAs dCooke et alj 120061) . (2) the identi fication of a 
numbe r of high-z DLAs associated with LB Gs ( M0ll er et"aT] 
2002a|b1:IChen et alJl2009l:lFvnbo et al.ll2^lFumagaliretai1 
201Qf lCookeetal]l2010h . (3) the occurrence of Ly-a emis- 
sion observed in the ce nter of DLA troughs (Mail er et alJ 
12004 ICooke et al.ll2010l) . and (4) the appearance of DLAs 
in the spectra of rare lensed LBGs at high redshift, where 
the DLA and LBG have similar redshifts dPettini et al.|[2 002b; 
iCabanac et alJl2008l ; lDessauges-Zavadskv et al.ll2010l) . 

In fact, recent results are consistent with the idea that gas in 
spatial ly extended DLAs encompasses compact LBGs. lErbl 
(2008) demonstrated that LBGs are rapidly running out of 
"fuel" for star formation, and cold (T ^100 K) gas in DLAs 
is a natural fuel source. This finding is supported by the 
measurement s of th e SFR and gas densities of LBGs by 
iTacconi et aD d2010|) . We note that if DLAs are the fuel 
source for the LBGs, the DLAs in turn would have to be 
replenished since the comoving density of DLAs at its peak 
(z ~ 3.5) is about 1/3 the curr ent cosmic mass density of 
stars (Prochaska & Wolfe 2009). Presumably, they are re- 
plenish ed through accretion of warm (T ~ 10 4 K) ioniz ed 
flowj 5 ! dDekel & Birnboimll2006L 12001; iDekel et al]|2009allbl; 
iBauermeister et aLll2010l) . 

While LBGs embedded in spatially extended DLA gas help 
resolve some properties of DLAs, it is problematic whether 
metal-enriched outflows from LBGs can supply the required 
metals seen in DLAs. Nor is it clear whether such outflows 
can generate turbulent kinetic energy at rates sufficient to bal- 
ance dissipative losses arising from cloud collisions implied 
by the multi-component velocity structure of DLAfl We note 

5 Often referred to as "cold" flows, but we call them warm flows since cold 
refers to T ~100K gas in this paper. 

The cloud crossing time is t cmss = H/v and the cloud collision time is 



that lFumagalli et all d201 ll) show that that the filamentary gas 
structures in the cold mode accretion scenario t hat provide 
galax ies with fresh fuel (e.g., De kel et aT]|2009at iKeres et ail 
2009) are not sufficiently dense to produce DLA absorption, 
nor do these filaments have a large en ough area covering frac- 
tion dFaucher-Giguere & K eres 201 1). The only location with 
sufficient covering fraction and high enough densities for the 
gas to become self shielded is in the vicinity around galaxies 
(2-10kpc). 

To address these issues, we adopt the working hypothesis 
that in situ star formation occurs in the presence of atomic- 
dominated gas in the outskirts of LBGs, similar to the outer 
disks of local galaxies. Since most of the atomic-dominated 
gas at high redshift is in DLAs, we assume t hat this is DLA 
gas. W e emphasize that while the past results dWolfe & Chenl 
2006) set sensitive upper limits on in situ star formation in 
DLAs without compact star-forming regions like LBGs, no 
such limits exist for DLAs containing such objects. 

For these reasons we search for spatially extended star for- 
mation associated with LBGs at z ~ 3 by looking for regions 
of low surface brightness (LSB) emission surrounding the 
LBG cores. We are searching for in situ star formation on 
scales up to ~ 10 kpc, where the detection of faint emission 
would indicate the presence of spatially extended star forma- 
tion. It would also uncover a mode of star formation hith- 
erto unknown at high z, and could help solve the dilemmas 
cited above. We test the hypotheses that (1) the extended 
star formation is fueled by atomic-dominated gas as probed 
by the DLAs, and (2) star formation occurs at the KS rate. 
That is, we consider whether star formation occurs in the out- 
skirts of LBGs, whether that star formation occurs in atomic- 
dominated gas, and at what SFR efficiency the stars form. 

This paper is organized as follows. In Section 2, we de- 
scribe the observations used, and in Section 3 we we identify a 
sample of compact LBGs at z ~ 3 to search for extended LSB 
emission around the compact cores of the LBGs. In Section 4, 
we describe the image stacking technique, measure a median 
radial profile of the extended LSB emission, and discuss pos- 
sible selection biases of the observations. In Section 5, based 
on the observed radial profile, we calculate the corresponding 
SFR surface density distribution and the sky covering fraction 
of the extended LSB emission. We then calculate the in situ 
SFR density and metal production in these extended LSB re- 
gions. In Section 6, we develop a theoretical framework to 
connect known DLA statistics to the observed surface density 
distribution of SFR in the outskirts of LBGs, and to obtain an 
empirical estimate of the star formation efficiency in distant 
galaxies. In Section 7, we discuss the impacts of our analy- 
sis in understanding the star formation relation in the distant 
universe, and in interpreting the observed low metal content 
of the DLA population. 

Throughout this paper, we adopt the AB magnitude system 
and an ((lM,^i\,h) = (0.3,0.7,0.7) cosmology. 

2. OBSERVATIONS 

We implement our search for spatially extended star forma- 
tion around z ~ 3 LBGs in the most sensitive high-resolution 

'coll ~ ' where H is the DLA scale height, v is the cloud velocity, n is the 
number density of clouds, and cr is the geometric cross section of the clouds. 
Therefore, the ratio of the cloud crossing time to the cloud collision time is 
■j 2 ^ ~ noH ~ r, where r is the optical depth. Consequently, if r > 1 as 
is observed, the clo uds would dissipate on a timescale s hort compared to the 
crossing time (e.g., McDonald & Miralda-Escude 1999). 
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Figure 1. Comparison of the number c ounts of z ~ 3 LBGs in 0.5 mag bins per square arcminute. The red circles are from Rafelski et al. 1 2009), the green crosses 
are from Sawicki & Thompson 12006) corrected for completeness, and the blue triangles are from Reddy & Steidel 12009) also corrected for completeness. The 
points from Reddv & Steidel (2009) are offset by 0.05 mag in order to avoid overlapping the Rafelski et al. (2009) points. The black line is the best-fit Schechter 
function I Schechter 1 976 ) from Reddv & Steidel HMk with a = -1.73 ± 0.13, Mt E (1700) = -20.97 ± 0. 14, and <p* = (1.71 ±0.53)x 10~ 3 Mpc~ 3 . 



images available: the V band image of the UDF taken with the 
Hubble Space Telescope (HST). This image is ideal because 
(1) of its high angular resolution (PSFFWHM = 0."09), (2) at 
z ~ 3 the V band fluxes correspond to rest-frame FUV fluxes, 
which are sensitive measures of SFRs, since short-lived mas- 
sive stars produce the observed UV photons, and (3) the la 
point source limit of V=30.5 implies high sensitivity. Since 
most of the LBGs in the UDF are too faint for spectroscopic 
identification, the u band is needed to find z ~ 3 LBGs via 
their flux decrement due to the Lyman limit through color se- 
lection and photometric redshifts. To this end, we acquired 
one of the most sensitive u ba nd images ever ob tained and 
identified 407 LBG s at z ~ 3 (Rafels ki et all 120091, see also 
iNonino et al.l d2009l) ). 

Throughout the paper we utilize the B, V, i', and z' band 
(F435W, F606W, F 775W, and F850LP, r espectively) obser- 
vations of the UDF dBeckwith et al.ll2006l) . obtained with the 
Wide Field Camera on the HST Advanced Camera for surveys 
(ACS: lFordetal.ll2002l) . These images cover 12.80 arcmin 2 , 
although we only use the c entral 1 1 .56 arcmin 2 which over- 
laps the u band image from (Rafelski et al. 2009). The u band 
image was obtained with the Keck I telescope and the blue 
channel of the Low-Resolution Imaging Spectrometer (LRIS; 
lOkeet al.|[l995t iMcCarthv et ail [19985 and has a ler depth 
of 30.7 mag arcsec" 2 and a limiting magnitude of 27.6 mag. 
The sample described below also makes use of the observa- 
tions taken with the NICM OS camera NIC3 in th e J and H 
bands (Fl 10W and F160W: lThompson et al.ll2005l) whenever 
the field of view (FOV) overlaps. 

3. SAMPLE SELECTION 

In order to search for spatially extended star formation asso- 
ciated with z ~ 3 LBGs, we require a sample of such galaxies 
to form a super-stack of LBG images that we describe here. In 
Section 3.1, we compare the number counts of the z ~ 3 LBGs 
in the UDF to those in the literature. Then in Section 3.2, we 
select a subsample that is appropriate for stacking in order 
to improve the signal-to-noise (S/N) in the LBG outskirts as 



described in Section 4. Lastly, in Section 3.3 we investigate 
possible selection biases of the observations. 

The samples in this paper are based on the z ~ 3 LBG sam- 
ple of iRafelski et al.1 (120091) . which contains 407 LBGs se- 
lected by using a combination of photometric redshifts and 
the u band drop out tech nique dSteidel & Hamilton! 119921 : 
ISteidel etal.|[l995l Il996allbl) . This selection of LBGs is en- 
abled by the extremely deep u band image described in Sec- 
tion 2, needed to reduce the traditional degeneracy of col- 
ors between z ~ 3 and z ~ 0.2 galaxies that c an yield in- 
correct redshifts at z ~ 3 w i thout th e u ba n d dEllisI 119971: 
| Fernandez-Soto et al] [19991: iBenftezI I2000t IRafelski et ail 
2009). Rafelsk ieTafl (120091) found that the resultant sample 
is likely to have a contamination fraction of only ~ 3%. Any 
such contamination will have a minimal effect on our results 
and is included in the uncer tainties (see Section 4.3). In ad- 
dition, IRafelski etall (12009b found that the LBG sample is 
complete to V ~ 27 mag, limited by the depth of the u band 
image. 

3.1. Number Counts 

In order to (1) be confident in our sample selection, (2) 
verify that we are probing the correct comoving volume, 
and (3) make completeness corrections in Sections 5 and 
6 (as described in Appendix A), we comp are the num- 
ber co unts of the z ~ 3 LBG selection from Rafelski et al. 
( 2009) to the comple t eness correc ted number coun t s from 
ISawicki & T hompson! (120061) and IReddv & Steidell (120091) 
in Fig ure Q] The uncertainties shown from IRafelski et alJ 
(2009) are only poisson and therefore h ave smaller error 
bars than those by Reddv & SteideJ (I2009I) . The points from 
IReddv & Steidell ( 12009b are off set by 0.05 mag for clarity, and 
their uncertainties include both poisson and field to field vari- 
ations. The numb er co unts and the best-fit Sch echter function 
dSchechteril 19761) from lReddv & Steidell d2009b are converted 
to number of L BGs per half-magn itude bin using a similar 
conversion as in Re ddv etail (2008). Specifically, we use the 
R band as a tracer of rest-frame 1700 emission. The apparent 
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Figure 2. Thumbnail images in the 
FigureF5] which corresponds to 18.5 



V band of the 48 z ~ 3 LBG subsample. The thumbnails are 2.4 arcsec on a side, the same size as the composite image in 
kpc at z ~ 3. 



measured R magnitude, Rab, is converted from the absolute 
magnitude using the relation: 



Rab = Mab ( 1 700) + 51og(<4/ 1 Opc) + 2 .51og( 1 + z) 



(1) 



where Mab (1700) is the absolute magnitude at the rest- frame 
1700, and is the luminosit y dista nce. We apply the k- 
correction from iRafelski et all (120091) of Rab - Vab= -0.15 
mag to get the V band mag. We use a comoving volume 
of 26436 Mpc 3 for the redshift interval 2.7< z <3.4 and 
an area of 1 1 .56 arcmin 2 for the number count conversion. 
We adopt a Schechter function with parameters found by 
iReddv & Steidel d2009l) of a = -1.73 ± 0.13, M AB (1700) = 
-20.97 ± 0. 14, and (f>* = ( 1 .7 1 ± 0.53) x 10" 3 Mpc" 3 . 

The resultant numb er counts agree nicel y, and show that the 
completeness of the Ra felski et al.l (120091) z ~ 3 LBG sam- 
ple matches the previous completeness limit found of V ~ 27 
magnitude. The agreement also suggests that the comoving 
volume for the redshift interval 2.7< z <3.4 is appropriate 
for the sample, which is important for the argument given in 
Section 5. More importantly, the agreement of the number 
counts allows us to use the best-fit Schechter function from 
IReddv &~S teidel (2009) to determine the expected number of 
LBGs at fainter magnitudes, providing the needed informa- 
tion to make completeness corrections in Appendix A. We 
note that this luminosity function is valid to R ~ 26.5, after 
which the extrapolation to fainter magnitudes could be a po- 
tential source of error, and we address this below. 

3.2. Catalogs of LBGs and Stars 

We use two catalogs of z ~ 3 LBGs in this paper. The 
first is the ful l samp le of 407 z ~ 3 LBGs as described in 
IRafelski et al.l (I2009I) . The sample redshi ft distribution is 
shown in Figure 12 of IRafelski etall (120091) and has a mean 
photometric redshift of 3.0 ± 0.3. The second (hereafter re- 
ferred to as "subset sample") is a sample of z ~ 3 LBGs se- 
lected to create a composite image to improve the signal-to- 



noise of the surface brightness profile described below. These 
LBGs are selected to be compact, symmetric, a nd isolated, 
similar to the selection done at higher redshift by Hathi et al. 
(2008). The LBGs are selected to be compact and symmetric 
to aid in stacking LBGs of similar morphology and physical 
characteristics such that the bright central regions of the LBGs 
overlap. They were also selected to be isolated from nearby 
neighbors to avoid coincidental object overlap and dynami- 
cally disturbed objects. 

To select objects that are compact, symmetric, and iso- 
lated, we measure morphologic al parameters in the V 
band image using SExtractor ( Ber tin & Arnoutsl 119961) . 
We experimented with d ifferent morpholog ical parame- 
ters, such as asymmetry ( Schad e~et al.l 119951) . concentra- 
tion dAbraham et al.ll 1994 1 19961), Gini c oefficient dLotz et all 



12004 ). and dumpiness dConsel ice 2003). However, we found 
that the best sample was selected based on the FWHM for 
compactness and ell ipticity e = (1 - b/a ) for symmetry, sim- 
ilar to the criteria in Hat hi et al.l (|2008). Specifically, for the 
subset sample, we require that FWHM < 0."25 and e <. 0.25 , 
a slightly more conservative selection than lHathi et al.l |2008). 
Lastly, to select isolated objects, we require that there are no 
other objects within l."4 brighter than 29th mag. These re- 
quirements yield a sample of 48 LBGs, representing ~ 12% 
of the z ~ 3 LBG sample, whose properties are compared in 
Section 3.3. We show this sample as thumbnails that are 2.4 
arcsec on a side, which corresponds to 18.5 kpc at z ~ 3, in 
Figure [2] and a table with relevant information about the 48 
LBGs in TableQ] Inf ormation on th e full sample of 407 LBGs 
is available in Rafels ki et al.l (12009). 

In addition to the two samples of LBGs, we also need 
a sample of stars for an accurate measurement of the point 
spread functi on (PSF) of the UD F images. We obtain this star 
sample from Pirzkal et al.l d2005l) . using only those stars with 
confirmed grism spectra, which mostly consist of M dwarfs. 
We exclude the stars that are saturated, leaving 15 stars in the 
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Table 1 

Properties of LBGs Included in Image Stack 



ID 11 


b 

£phot 


V 


u-V 


B-V 


V-z' 


FWHM 


Ellipticity 






(mag) 


(mag) 


(mag) 


(mag) 


(arcsec) 




84 


Til TU.'HJ 

J - 11 -0.40 


26.56±0.02 


2.40±0.35 


0.70±0.05 


0.03±0.04 


0.24 


0.18 


862 




27.16±0.02 


1.87±0.36 


0.74±0.05 


-0.44±0.05 


0.12 


0.14 


906 




27.52±0.02 


1.50±0.36 


-0.10±0.04 


-0.14±0.05 


0.10 


0.04 


1217 


2 7 =+0.37 


26.53±0.01 


2.28±0.30 


0.34±0.02 


0.01 ±0.02 


0.14 


0.23 


1273 


3.03^1 


26.24±0.01 


2.79±0.37 


0.58±0.03 


0.02±0.02 


0.12 


0.12 


1414 




27.19±0.02 


1.82±0.37 


0.49±0.06 


-0.03±0.05 


0.10 


0.13 


1738 


2 67 +0.36 

3 44+8:2f 

J - 4i *-0.43 


26.27±0.01 


1.05±O.O8 


0.03±0.02 


-0.30±0.03 


0.16 


0.15 


1753 


27.47±0.03 


1.02±0.35 


1.39±0.14 


0.58±0.04 


0.18 


0.12 


2581 


•3 4(1+0-43 
J -^ u -0.43 


26.92±0.02 


2.05±0.35 


1.04±0.07 


0.43±0.03 


0.18 


0.24 


2595 


2 Q7+0.39 


27.37±0.02 


1.60±0.35 


0.22±0.04 


-0.17±0.05 


0.13 


0.05 



Note. — V magnitudes are total AB magnitudes, and colors are isophotal colors, u band photometry is from 
Rafelski et al. ( 2009) and the rest are from Coe et al. 1 2006). Non-detections in the a band are given 3a limiting 
magnitudes. This table is available in its entirety in a machine-readable form in the online journal. A portion is 
shown here for guidance regarding its form and content. 

11 ID numbers from Rafelski et al. 1 2009) which match those by Coe et al. i 2006). 
Bayesian Photometric Redshift (BPZ) and uncertainty from 95% confidence interval from Rafelski et al. ( 2009). 



V band image within our FOV with V = 25.7 ± 1 .3, more than 
adequate to measure the PSF. We note that a comparison of 
the star PSF with the LBGs shows that they are all resolved in 
the high resolution ACS images. 

3.3. Comparison of the Subset and Full LBG Samples 

We investigate whether the subset sample of 48 LBGs is 
drawn from the same parent population of the full sample 
of 407 LBGs by comparing the magnitude, color, and red- 
shifts of two the samples. First, we find little variation in the 
magnitude distributions of the two samples, with a difference 
in the mean of ~ 0.3 mag. The subset sample is somewhat 
fainter, with the full sample having an average AB magnitude 
of V = 26.4 ± 0.9 and the subset sample with V = 26.7 ± 0.6. 
That is, there is a minor systematic selection of fainter LBGs 
in the subset sample, although this difference is not signif- 
icant. The similar magnitude distribution of the rest-frame 
FUV flux suggests that the SFR of the two samples is similar. 

Second, we compare the mean colors of the of the two sam- 
ples and find the two samples have the same colors. We test 
this both for the distribution and for stacks of the LBGs. First, 
the mean of the distribution of the subset sample yields col- 
ors of B-V = 0.5±0.4, V-i' = 0.0±0.2, and i'-z' =0.0±0.1, 
while the full sample has colors of B-V = 0.6±0.4, V-i' = 
0.1±0.2, and i'-z' =0.0±0.1. Second, the color of the stacked 
subset sample based on aperture photometry has colors of 
B-V = 0.2±0.2, V-i' = 0.1 ±0.2, and i'-z' =0.1 ±0.2, while 
the full sample stack has colors of B-V = 0.3±0.1, V-i' = 
0.1 ±0.1, and i'-z' =0.2±0.1. These colors are not signifi- 
cantly different based on both the distribution and the stacked 
photometry uncertainties. The similar distribution of colors 
suggests that the two samples are made of the same stellar 
populations and that their star formation histories (SFHs) are 
similar. 

Lastly, the two samples have very similar redshift distribu- 
tions, with the same mean redshift of 3.0 ± 0.3. We therefore 
conclude that the subset sample and full sample of LBGs are 
equivalent in magnitude, color, and redshift, and therefore are 
probably drawn from the same parent population of LBGs that 
have similar SFRs, stellar populations, and SFHs. Hence, we 
are relatively confident that the results determined below for 
the subset sample of LBGs are applicable to the full sample. 



For the sake of completeness, we also consider stacking the 
full stack of LBGs in Appendix B. We note that stacking the 
full sample introduces contamination into the stack, such as 
nearby galaxies. In addition, bright parts of morphologically 
different galaxies contribute to the faint parts of other galax- 
ies. We therefore do not use this as our primary stack, and 
take the full stack as an upper limit to the emission from the 
full sample of LBGs. 

4. ANALYSIS OF UDF IMAGES 

The V band image of the UDF is the most sensitive high 
resolution image covering the rest-frame FUV at z ~ 3 avail- 
able. However, even this image does not reach the desired 
sensitivity to search for spatially extended star formation on 
scales up to ^10 kpc (as shown below). We wish to in- 
crease the S/N high enough to probe down to low values of 
the SFR surface density (Esfr)- Image stacking methods can 
be used to study the average properties of well defined sam- 
ples i n which individual objects do not have the necessary S/N 
(e.g.. iPascarelle et al.|[l996t IZibetti etafl 12004 120051 120071: 
Hath i et al.ll2008l) . 

We therefore create super-stacks of the LBG images in Sec- 
tion 4.1 and investigate how the sky-subtraction uncertainty 
affects those stacks in Section 4.2. Using the super-stack of 
the subset sample of LBGs described in Section 3.2, we de- 
termine the radial surface brightness profile of z ~ 3 LBGs 
in Section 4.3, which will be used for much of the analysis 
throughout the paper. Lastly, in Section 4.4, we investigate 
the effects of the Lya line on the stacked image. 

4. 1 . Image Stacks 

We create stacked composite images for the LBG subset 
sample, full LBG sample, and stars sample using custom IDL 
code. For e ach object, we fit a two-dimensional Gaussian us- 
ing MPFIT dMarkwardj 120091) to determine a robust center. 
We then shift each object to be centered with sub-pixel res- 
olution, interpolating with a damped sine function. We then 
create thumbnail images for each object and combine all the 
objects by taking the median of all the thumbnails, yielding a 
robust stacked image which is not sensitive to outliers. While 
some of the individual thumbnails may have faint emission 
regions below the level of the isolation criteria, they do not 
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Figure 3. Composite image in the V band of the 48 z ~ 3 LBG subsample. 
Each stamp is 2.4 arcsec on a side, which corresponds to 18.5 kpc at z ~ 3. 



contribute to the median because such emission would need 
to occur at the same pixels for a significant number objects 
to affect the median. Therefore, independent of the origin of 
any such faint regions, they do not affect their median, and 
therefore do not affect the final stack. While we are most in- 
terested in the V band, we also carry out this procedure for the 
B, V, i', and z' bands and the stars sample. Figure[3]shows the 
stacked V band image, as this corresponds to the rest-frame 
FUV luminosity at z ~ 3 which is a sensitive measure of the 
SFR. The image is 2.4 arcsec on a side, which corresponds to 
18.5 kpc at z~ 3. 

We showed in Section 3.3 that the subset sample is rep- 
resentative of the full LBG sample and here investigate any 
possible variations in the radial surface brightness profile with 
magnitude and FWHM within the subsample itself. We check 
if these parameters affect the stack by creating three indepen- 
dent stacks of different brightnesses or FWHM. In the case of 
magnitude, we create one stack of the brightest 16 LBGs, a 
second stack of the next 16 LBGs, and a third stack with the 
faintest 16 LBGs, and repeat for FWHM. We find the profiles 
to be very similar and find that the magnitude and FWHM 
range does not affect our stack. We also investigate the change 
in the surface brightness profile color. We find that the vari- 
ations of the LBG composite images for the B, V, i', and z' 
bands across radius are small compared to their uncertainties, 
and no clear change in color is obvious at any radius. 

We also test the difference between taking the median and 
the mean of the images in our stack. The profiles of the mean 
stack are slightly higher at the bright end, but are very similar 
to the median stack starting at ~0."3. We chose to go with the 
median as it is more robust to possible contamination in the 
outskirts. In this way, we are not sensitive to contamination if 
it only occurs in a small subset of our sample. 

In creating a stack of the star data to measure the PSF, we 
scale each star to the peak of the fitted two-dimensional Gaus- 
sian before stacking. For the PSF, we only care about the 
shape of the PSF, and not the actual value of the flux. Given 
the wide distribution of magnitudes of the stars, this scaling 
improves the PSF determination. We do not scale the LBGs 
before stacking, as we care about the actual measured flux in 
the outskirts of the LBGs. We find that scaling does not have 
a significant effect on the LBG stack due to the small range 
in brightnesses of the selected LBGs, with the radial surface 
brightness of the two stacks being consistent within the un- 



certainties, and this would therefore not alter our results. 

4.2. Sky -subtraction Uncertainty 

We hope to accurately characterize the sky background 
and the uncertainties due to the subtraction of this sky back- 
ground. The sky-subtra ction uncertainty o f the UDF was 
carefully investigated by Hat hi et ail ((2008), and we follow 
their prescription for de termining the la sky-subtraction er- 
ror. lHathi et al.1 d2008l) found that it was more reliable to 
characterize the sky locally rather than globally for the en- 
tire image. For this reason, we redetermine the local sky 
background for the 407 z ~ 3 LBGs in our sample. First, 
we measure the sky background in each of the thumbnail 
images of the full LBG sample using IDL and the proce- 
dure MMM.pr cQ, adapted fro m the DAOPHOT routine by the 
same name dStetsonl 119871) . This procedure iteratively de- 
termines the background, removing low probability outliers 
each time, until the sky background is determined. We then 
find the la sky value by fitting a Gaussian to the distribu- 
tion of sky values. For the V band, this gives us a value for 
c s ky,ran of 3.52 x 10~ 5 electrons s" 1 , very similar to that found 
bv lHathi et ail d2008l) of 3.5 5 x 10" s electrons s " 1 . Using this 
number and the formalism in lHathi et al.l (120081) . we find a la 
sky-subtraction error of 30.01 mag arcsec" 2 . 

If the error is random, the uncertainty of the sky-subtraction 
will decrease as we stack more images together. In fact, for 
a median stack in the Poisson limit, the la uncertainty is 
1.25/vW, where N is the number of images. We test this 
relation specifically for the UDF data as the error may not 
be completely random. We stack 48 blank thumbnails and 
compare the standard deviation of stacked pixels to the me- 
dian of the standard deviations of each individual thumbnail 
and find the above relation holds to 99% accuracy. We are 
therefore confident that the sky-subtraction noise decreases 
with stacking as expected, yielding larger S/N values. For the 
stack of 48 LBGs from our subset sample, we find a la sky- 
subtraction error of 31.87 mag arcsec" 2 . We use both of these 
sky-subtraction errors as our sky limits below in Section 4.3. 

4.3. Radial Surface Brightness Profile 

We extract a radial surface brightness profile from the 
super-stack of LBG images in Figure [4] using custom IDL 
code which yields identical results to the IRAF3 procedure 
ELLIPSE. We use circular aperture rings with radial widths 
of 1.5 pixels in such a way that they do not overlap to avoid 
correlated data points while still finely probing the profile. 
We use custom code in order to facilitate using the bootstrap 
error analysis method to determine the uncertainties. Since 
we are stacking different objects with different characteris- 
tics, the uncertainty for the brighter regions will be dominated 
by the sample variance, and can be determined with the boot- 
strap method. Specifically, we bootstrap to get the uncertainty 
of the median of the LBGs by replacing a random fraction 
(1/esi 37%) of the 48 LBG thumbnails with randomly du- 
plicated thumbnails from the subset sample. We repeat this 
1000 times to get a sample of composite images using the 
method described in Section 4.1, each with its own radial sur- 
face brightness profile. The resultant uncertainty is then the 

7 Part of the IDL Astronomy User's Library 

8 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 
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Figure 4. Extracted surface brightness profile from the stacked image, where 
the black points are the 48 LBG composite profile and the blue line is the 
measured point-spread function from stars. The dotted red line is the best- 
fit Sersic profile convolved with the PSF. The dashed line is the lcr sky- 
subtraction error for one thumbnail image, while the dash-dotted line is the 
lcr sky-subtraction error for the composite stack of 48 LBGs, as described in 
Section 4.2. 

standard deviation of all the surface brightness magnitudes 
at each radius. This method is conservative, but includes all 
the uncertainties associated with the variance of the sample. 
Moreover, it helps takes into account possible errors intro- 
duced by possible contamination of our LBG sample by other 
objects, although as mentioned in Section 3, we believe this 
contamination fraction to be small. 

The solid blue line in Figure|4]represents the measured ACS 
V band PSF determined from a stacked image of the 15 stars 
described in Section 3.2. The central surface brightness of the 
stars is scaled to match the LBGs. We note that the S/N of the 
star stack is higher than that of the LBG stack, even though it 
has a smaller number of objects in the stack, because the stars 
are brighter than the LBGs. Therefore the PSF is well deter- 
mined, and we expect the uncertainties are dominated by the 
lower S/N LBG stack. The PSF declines more rapidly than 
the radial surface brightness profile at all radii, which shows 
that the LBGs are clearly resolved and that the median sur- 
face brightness profile of the LBGs is extended. The dashed 
line is the lcr sky-subtraction error for one thumbnail image, 
while the dash-dotted line is the lcr sky-subtraction error for 
the composite stack of 48 LBGs, as described in Section 4.2. 

4.3.1. Sersic Profile Fit 

The dotted red line in Figure |4] shows the best fit Ser- 
sic profile convolved with the ACS V band PSF. We fit for 
the best Sersic model to the composite LBG image by using 
the Levenberg-Marquardt least-squares minimization, with 
the x 2 calculated on a pixel by pixel basis for each pos- 
sible model. The best fit values are n = 1.9 ±0.04 and 
R e = 0.074 ±0.001 arcsec, where n is the Sersic index and 
R e is the effective radius, which includes half the light of 
the LBGs. We use the dimensionless scale factor b(n) from 
ICiotti & Bertinl (119991) such that R e is the half-light radius. 
The small R e value is likely due to our pre-selection of LBGs 
to be compact (see Section 3.2). While the fit has a good re- 
duced chi-square x 1 l v °f 1-27, it is not a very good fit to the 
data in the outer regions. The fit is dominated by the inner 
part of the profile with radii less than ~0."4, since the uncer- 
tainties are significantly smaller in that region. For radii larger 
than ~0."4, the profile deviates from the inner best-fit profile. 



This deviation is real, being above the PSF and the la sky- 
subtraction error. This is similar to what lHathi et ail (120081) 
found when stacking LBGs at z ~ 4 — 6. The main constraint 
we have from this is that it suggests that the profile is similar 
to an exponential disk type profile. If we fit an exponential, it 
yields a worse fit with a x l v of 1-4, and R e = 0.068 ±0.001 
arcsec. A bulge-disk model yields a slightly better fit in the 
outer regions, but does not improve the overall x 1 l v - We note 
that we do not use the fit in the analysis below. 

4.4. Effects of the Lya Line on the Image Stack 

We investigate possible contamination from Lya emission 
on the image stack to ensure that the radial surface brightness 
profile is unaffected by Lya emission from the LBGs. First, 
we note that Lya only enters our V band filter for about half 
our sample due to the redshift range sampled. Second, we find 
that due to the large width of the V band filter, the Lya line 
would have a very small effect. This is determined by tak- 
ing the stacked LBG spectrum from Sha plev et alj d2003l) . and 
comparing the flux in the V band filter with and without the 
Lya line. We find that the inclusion of the Lya line yields an 
increase of 0.02 mag, a very small effect compared to our un- 
certainties. While our sample of LBGs is signi ficantly fainter 
than th e LBGs in the stacked spectrum from Shaple yet al.1 
(2003), we expect the average strength of the Lya line to be 
similar due to the low escape f raction of ionizing radiation 
from LBGs (IShaplev et al.ll2006l) . 

Moreover, we compare the radial surface brightness profile 
of the V band stack and an equivalent stack in the i' band. The 
Lya line does not enter the i' band filter throughout our red- 
shift range, making the i' band an excellent test to check if the 
radial surface brightness profile is affected by the Lya line. 
We find that the V and i' bands have the same radial surface 
brightness profile within their uncertainties, with no system- 
atic shifts. This suggests that the Lya line has little to no 
effect, even if the Lya line were more extended than the con- 
tinuum. 

5. DIRECT INFERENCES FROM THE DATA 

The surface brightness profile of the stacked image indi- 
cates the presence of spatially extended star formation around 
LBGs. In this section, we describe how to connect this emis- 
sion, which corresponds to the rest-frame FUV radiation in- 
tensity, to the SFR surface density. We then compare the 
covering fraction of this emission to that of the gas respon- 
sible for the star formation in Section 5.2. Specifically, in 
order to determine what types of gas can be responsible for 
the observed emission, we compare its covering fraction to 
that of atomic-dominated gas in Section 5.2.1 and to that of 
molecular-dominated gas in Section 5.2.2. Then in Section 
5.3 and 5.4, we calculate the in situ SFR, p\, and metal pro- 
duction based on the integrated flux measured in the outskirts 
of LBGs. 

5.1. Connecting the Observed Intensity to the SFR Surface 

Density 

In order to investigate how the measured star formation re- 
lates to the underlying gas, we require a relation between star 
formation and gas properties. Star formation occurs in the 
presence of cold atomic and/or molecular gas, according to 
the KS relation given by 

£sFR = *x0g^. (2) 



8 



Rafelski, Wolfe, & Chen 



This relation holds for nearby disk galaxies in which E gas is 
the mass surface density perpendicular to the plane of the disk, 
E c = lM Q pc" 2 , K = KK erm =(25±0.5)x \Q- 4 Mq yr" [ kpc" 2 , 
and /3=1.4±0.15 (iKennicuttl Il998allbb - There has been 
much recent work on improving both our understanding of 
this relation, and measuring the values of K and 8 (e.g. , 
Lerov et al.ll2008t iBigiel et al.ll2008t iKrumholz et aljr2 009b: 



Gnedi n & Kravt sov 20101: iGenzelet all 12010; BigielejtalJ 
201 Obi) . When only considering molecula r gas, the relation- 
ship has a flatter slope of /?=0 .96±0.07 dBigiel et al.l 120081) 
or 8 ~ 1.1 (Wong & Blitz 2002). We use the original values 
from ( Kennicutt 1 1 998al) when considering the total gas den- 
sity to simplify comparisons with other work, and 8=1.0 and 
#=#Biegei=(8.7±1.5)X 10" 4 M© yr" 1 kpc" 2 when considering 
only molecular gas. We note that t he K value given he re used 
for molecular gas is modified from Bigi el et all d2008l) to use 
the same S c = 1 M Q pc" 2 value as above. 

Rewriting the KS relation in terms of the column density of 
the gas, we get 



N 



(3) 



where the scale factor iV c =1.25x 10 20 cm" 2 (Ke nnicuttl 
Il998allbh and N is the hydrogen column densityO We note 
that this is only valid above the critical column density, which 
is usually associated with the threshold condition for Toomre 
instability. For H I gas in local galaxies, it is observed to range 
between 5 x 10 20 cm" 2 and 2x 10 21 cm" 2 (lKennicuttlll998bl) . 

In order to connect Esfr to the observations, we require 
a relation between observed intensity, corresponding to rest- 
frame FUV emission, and ob served column density, N. Fol- 
lowing lWolfe & ChenI (12006. equation (3)), we find that for a 
fixed value of N, the intensity averaged over all disk inclina- 
tion angles is given by 



robsx 



SFR 



4tt(1+z) 3 /3 ' 



(4) 



where z is the redshift, and C is the conversion fac- 
tor from SFR to FUV (A ~ 1500) radiatio n, with C = 
8.4 x 10~ 16 erg cm" 2 s" 1 \izr l (M P) yr" 1 kpc" 2 )" 1 dMadau et alj 



1998t IKennicuttl Il998bl) . We use the same value of C as 
Wolfe & ChenI d2006l) corresponding to a Salpeter IMF. The 
result in Equation assumes that the star formation occurs 
in disks inclined on the plane of the sky by randomly selected 
inclination angles an d averages over all possible angles (see 
IWolfe & Chenfl2006l) . 

5.2. Covering Fraction ofLBGs Compared to the 
Underlying Gas 

The covering fraction of observed star formation should be 
consistent with that of its underlying gas. We therefore inves- 
tigate whether the covering fraction of the outer parts of LBGs 
is consistent with the covering fraction of atomic-dominated 
gas in Section 5.2.1 and molecular-dominated gas in Section 
5.2.2. This consistency check yields insights into the nature of 
the observed star formation and validates the hypothesis that 

9 The reader should be aware that when referring to nearby galaxies, N cor- 
responds to N± , the H I column density perpendicular to the disk, but when 
writing about our observations, we are referring to observed column densities 
N, where we implicitly include the inclination angles in our definitions. 



the outskirts of LBGs consist of atomic-dominated gas, which 
is used in subsequent sections of the paper. 

We calculate the cumulative covering fraction, Ca, for gas 
columns greater than some column density, N, by integrating 
the H column-density distribution function f(Nn,X), where H 
is either H lor H2. Specifically, 



C A (N) = 



dX / dN H f(N H ,X) . 



(5) 



where f(Nn,X) is the observed column-density distribution 
function of the hydrogen gas, N max is the maximum column- 
density considered, and X is the absorption distance with dX 
being defined as 



dX 



H ° {l+z?dz, 



H(z) 



(6) 



where Hq is the Hubble constant and H(z) is the Hubble pa- 
rameter at redshift z. For X m \ n and X m . dx we use the same red- 
shift interval as in Section 3.1, namely 2.7 < z < 3.4 corre- 
sponding to 6.6 < X < 9.2. The covering fraction depends 
strongly on the column-density distribution function, which 
is different for atomic and molecular gas. Below we investi- 
gate the covering fraction for both cases. 

5.2.1. Covering Fraction of Atomic-dominated Gas 

There is strong evidence to support the association of LBGs 
and neutral atomic-dominated H I gas, i.e., DLAs (see Section 
1), and we therefore investigate whether the covering fraction 
of the outer parts ofLBGs is consistent with the covering frac- 
tion of DLAs. If the outer regions of LBGs truly consist of 
DLA gas, then the covering fractions as functions of the sur- 
face brightnesses should be consistent. In this subsection, we 
work under the hypothesis that the observed FUV emission in 
the outskirts of LBGs is from in situ star formation in atomic- 
dominated gas and compare the observed covering fraction to 
that of the gas distribution. 

In order to calculate the covering fraction using Equation 
(O, we require f{Nu,X) for atomic-dominated gas. The ob- 
served H I column-density distribution function, f(Nm,X), is 
obtained by a double power-law fit to the Sloan Digital Sky 
Survey data: 



-24 2 

cm , a=ai,= 



(7) 



-2.00±0.05 for 



where fc 3 =(1.12±0.05)x 10 

N < JV(H a nd g=a 4 =-3.0 for N > N , where N = 3.54^ : 
10 21 cm" 2 dProchaska et al J 120051: iProchaska & WoIH 20b9). 
The value of as, used is different than measured in 
dProchaska & W olfe 2009), to remain consistent with our for- 
mulation of randomly oriented disks in Section 6, and is pre- 
dicted to be -3.0, and we use this value for the covering frac- 
tion to be consist ent. Although this value o f ct\ is different 
than the value in Prochaska & Wolfe (2009), the uncertain- 
ties are quite large due to low numbers of very high column 
density DLAs, and it is quite similar to the value found by 
iNoterdaeme etall d2009b of cu=-3 .48. 

We note that INoterdaeme eTall d2009l) find slightly differ - 
ent values for £3 and 0:3 than IProchaska & Wolfd d2009l) . 

10 We follow lWoIfe & Cheri f2WB) who equated N with N d . the break in 
the double power-law expression for fXNm.X). 
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Figure 5. Cumulative covering fraction of DLA systems with columns 
greater than some column density N (and therefore surface brightness), com- 
pared to the covering fraction of the observed z ~ 3 LBG outskirts greater 
than some surface brightness. The blue dotted line is the DLA covering frac- 
tion for DLAs with K = Ajcerm, me red dash-dotted line is for K = 0. 1 X A'Kcnn, 
and the red triple-dot-dashed line is for K = 0.02 X K^ am . The solid black line 
is the covering fraction for the LBG outskirts, while the gold long-dashed line 
is the same corrected for completeness. The gray filled region represents the 
transition region mentioned in Section 5.2.2. The top jc-axis coordinates are 
the column densities corresponding to the surface brightnesses using an effi- 
ciency of K = 0. 1 X XKenn ■ 



with £ 3 =8.1 x 10 cm 2 and a=«3=-1.60 for N < No, cor- 
responding to a flatter s lope. We use the values from 
Prochaska & Wolfe (2009) and describe how the differing val- 
ues affect our results below. Also, although the normalization 
of f(Nni,X) varies with redshift, the variations for our red- 
shift in terval are not large and do n ot strongly affect Ca- Us- 
ing the iProchaska & Wolfe! d2009l) values for f(N Hh X), and 
N mdx = 10 22 cm" 2 , we calculate the covering fraction using 
Equation (0 and the corresponding expected /ly from Equa- 
tions © and ©, where fi v = -2.51og(/° o bs ) -Z, and Z is the 
AB magnitude zero point of 26.486. 

We compare the cumulative covering fraction of DLAs to 
the observed covering fraction of the outer regions of the 
LBGs in Figure [5] The blue dotted line is the DLA covering 
fraction for DLAs forming stars according to the KS relation, 
with K = ATxenn- The red lines are for less efficient star forma- 
tion where the dash-dotted line represents K = 0. 1 x ^Kenn and 
the triple-dott-dashed line represents K = 0.02 x ^Kenn- The 
black line is the covering fraction for the observed emission 
in the outskirts of LBGs in the UDF, which is just the area 
covered by the outskirts of LBGs up to /iy divided by the to- 
tal area probed, 11.56 arcmin 2 . The covered area is obtained 
from the radial surface brightness profile and the 407 observed 
LBGs in this area. The observed covering fraction depends on 
the depth of the images and therefore requires a completeness 
correction for faint objects that are missed: we discuss this 
completeness correction in Appendix A. The completeness 
corrected covering fraction for LBGs is the dashed gold long- 
dashed line, and is the curve that we compare to the DLA lines 
below. We ignore the data at a radius > 0.8 arcsec (~ 6kpc) to 
only include data with S/N>3, however, we expect a contin- 
uation of the observed trends. We note that we are only sam- 
pling the top end of the DLA distribution function, and there- 
fore the covering fraction shown is a small fraction (about 
a tenth) of the total covering fraction of DLAs. When con- 
sidering all DLAs, we cover about one-thi rd of the sky, i.e., 
log(C/i) ~ -0.5. We note that if we use the iNoterdaeme et all 



(2009) values for f(N m ,X), then the DLA lines move up and 
to the left (i.e., cover more of the sky). 

Under the hypothesis that the outskirts of LBGs are com- 
prised of atomic-dominated gas (i.e., DLAs) which is respon- 
sible for the observed emission, then we expect the covering 
fraction of the DLA gas to be equal to the covering fraction of 
the outskirts of LBGs. The blue dotted line for DLAs fol- 
lowing the KS relation would therefore only be consistent 
with the covering fraction of the outskirts of LBGs if much of 
the DLA gas is not surrounding LB Gs. This possibility was 
constrained by (Wolfe & Chen 2006), who found that DLAs 
would need to be forming at significantly lower SFR efficien- 
cies if this was the case. On the other hand, we find that the 
covering fraction of the outskirts of LBGs is consistent with 
DLAs having SFR efficiencies of K > 0.1 x K^ enn , at which 
point the covering fraction is roughly equal. This covering- 
factor analysis provides evidence that if the outskirts of LBGs 
are comprised of DLA gas, then the SFR efficiency of this 
atomic-dominated gas at z ~ 3 is about 10% of the efficiency 
for local galaxies. 

Moreover, the cumulative covering fraction shows that 
there is sufficient DLA gas available to be responsible for the 
emission in the outskirts of LBGs for SFR efficiencies > 10%. 
We investigate this lower SFR efficiency further in Section 6, 
where we find an efficiency closer to 5%, which is only a fac- 
tor of ~2 different than the efficiency determined from the 
covering fraction. We note that systematic uncertainties due 
to assumptions throughout both quantities could easily be off 
by a factor of two, and so the general agreement of the cover- 
ing fraction at low SFR efficiencies is reassuring, and we are 
not concerned about the minor disagreement. 

5.2.2. Covering Fraction of Molecular Gas 

In the previous subsection, we worked under the hypothe- 
sis that the observed FUV emission in the outskirts of LBGs 
is from in situ star formation in atomic-dominated gas. How- 
ever, it is possible that this star formation occurs in molecular- 
dominated gas. We consider this scenario here, and compare 
the covering fraction of molecular-dominated gas, where the 
majority of the hydrogen gas is molecular, to our observa- 
tions. 

In order to calculate the covering fraction using Equation 
(0, we require f(Nn,X) for molecular-dominated gas. We 
use the observed molecular column-density di stribution func- 
tion, /(Nh 2 ), from lZwaan & Procha ska (2006), w ho obtaine d 
a lognormal fit to the BIMA SONG data dHelfer et alJ l2003). 



f(N H2 )=r exp 



logN-pi. 



/2 



(8) 



where /i = 20.6, a = 0.65, and the no rmalization /* is 1 . 1 x 
10 cm (IZwaan & Proc haska 20060 We use the molecu- 
lar version of the KS relation discussed in Section 5, where 
13=1.0 and ^=^ B ie g ei=8.7x 10" 4 M yr _1 kpc -2 , and we let 
N m - dx = 10 24 c m" 2 , the largest observed v alue for the f(Nu,) 
function used (Zwaan & Prochaska 2006). However, f(Nn,) 
is not observationally determined at z ~ 3, and likely evolves 
over time, and we investigate such a possibility below. 

The evolution of f(Nn,) would either be due to a change 
in the normalization or a change in the shape. The shape of 



1 1 We note that Zwaan & Prochaska 1 2006) have a typographical error, 
switching /i and u. 
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Figure 6. Cumulative covering fraction of gas with columns greater than 
some column density N (and therefore surface brightness). This figure is sim- 
ilar to Figure [5] but makes comparisons to the covering fractions of molecu- 
lar gas rather than atomic-dominated gas. The solid black line is the covering 
fraction for LBGs starting from the center of the LBG core, and the gold long- 
dashed line is the same corrected for completeness. The purple short-dashed 
line is the covering fraction of molecular hydrogen with no evolution at z ~ 3, 
the pink triple dot-dashed line is the same with an evolution of four times /* 
(the normalization of the column-density distribution function of molecular 
gas) and the cyan dot-dashed line is the same with an evolution of 10 times 
f* . The gray lines continuing the purple, pink, and cyan lines are extrapola- 
tions of the data to lower column densities. The gray filled region represents 
the transition region mentioned in Section 5.2.2. The top jc-axis coordinates 
are the column densities corresponding to the surface brightnesses using an 
efficiency of K = Xrjiegel an d ft = 1. valid for molecular gas. 

the atomic gas column-density distributio n function, /(NgX 
has not evolved between z = and z = 3 (jZwaan et al.ll2005t 
iProchaska et alJl2005l:|Pro"chaska & Worfefl2009l) . but the nor- 
malization has increased by a factor of two. In the case of H2, 
we consider the instance in which only the normalization (/*) 
evolves, then we are looking for a change in J1h, at z = to 
fln 2 at z = 3. While theoretical models predict that fl^(z = 
3)/f2ei(z = 0) is ~ 4, their s imilar prediction for atomic ga s 
does not match observations (lObreschkow & R awlings 2009). 
Alternatively, we can determine an upper limit of the evo- 
lution of /* using the evolution of p\ for galaxies between 
z = and z = 3, assuming that the evolution in f(N^) is only 
due to the normalization, and there is no evoluti on in the KS 
law for molecular gas between z = and z = 3 dBouche et al.l 
2007; Dad di et alJlfoToHGenzel et al.ll201d see Section 6.3). 
Specifically, studies have found that p'Jz = 3) //c4(z = 0) ~ 10 
dSchiminovich et aT1l2005t iReddv et ail 120081) . and therefore 
/* changes at most by a factor of 10, if we assume that the 
contribution from atomic gas is small. This is used as an up- 
per limit to the evolution of /* , and we are not suggesting that 
this is the correct evolution. 

We plot the covering fractions of molecular gas at z = 
3 in Figure [6] We consider three cases for /(TVh,): (1) 
no evolution as a purple short-dashed line, (2) evolution 
with a factor of four increase i n /* based on the model by 
Obreschkow & Rawlings (2009) as a pink triple-dashed line, 
and (3) evolution with a factor of 10 increase in f* based on 
the observed evolution in p\ as a cyan dot-dashed line. The 
gray lines continuing these three lines are extrapolations of 
the data to lower column densities than observed. We note 
that the column densities on the top jc-axis of the plot now are 

for A' = A' B iegel- 

These covering fractions of molecular gas are compared to 
the LBG profile including the inner cores in Figure [6] The 



LBG covering fraction is now modified to include the LBG 
cores which are composed of molecular-dominated gas, and 
the new LBG covering fraction is plotted as a solid brown line. 
We again ignore the data at a radius > 0.8 arcsec (~ 6kpc) to 
only include data with S/N>3. This covering fraction is also 
corrected for completeness similar to Section 5.2.1, except 
this time we include the cores of the LBGs and make no dis- 
tinction between the outskirts and the inner parts of the LBGs, 
and is described in Appendix A. 1 . We plot this corrected cov- 
ering fraction as a gold long-dashed line. This correction is 
quite large, as even though the cores cover a smaller area than 
the outskirts, the cores of fainter missed LBGs contribute at 
every surface brightness as we go fainter. Since there are sig- 
nificantly more faint LBGs than bright LBGs, there are signif- 
icantly more LBG cores contributing to each surface bright- 
ness than there are LBGs with outskirts at those same surface 
brightnesses. This correction assumes that all the star forma- 
tion comes from molecular-dominated gas, and therefore all 
the cores of fainter LBGs are included. 

As in Section 5.2.1, we expect the covering fraction of the 
gas to be equal to the covering fraction of the LBGs forming 
out of that gas. In the case of purely molecular gas, the upper 
limit of the covering fraction (10 x /*) of the gas is only con- 
sistent to pv ~ 28.5. At pv > 28.5, the covering fraction of 
LBGs is larger than that of the molecular gas. In order to have 
a covering fraction larger than all the corrected LBG emis- 
sion, we would require an evolution of /* by a factor of more 
than 60, which is not very likely. 

We have just shown that the predicted covering factor for 
molecules is inconsistent with our results for py > 28.5. At 
the same time, our results for atomic -dominated hydrogen 
do not apply for py brighter than ~ 29, due to the atomic- 
dominated H I gas cutoff of Nm < 10 22 cm" 2 . Therefore, the 
surface brightness interval from 28.5 < pv < 29.2 is not si- 
multaneously consistent with the neutral atomic-dominated 
gas nor the molecular gas' covering fractions. This lower 
bound on pv is based on comparing the LBG data compared 
to 10 times /*, and for lower evolutions of /*, this begins 
at even brighter py. This result is reasonable if the LBG 
outskirts consist of atomic-dominated gas. In this scenario, 
there would need to be a transition region between atomic- 
dominated and molecular-dominated gas, where star forma- 
tion occurs in both phases. We note that in this hybrid re- 
gion, the corrected covering fraction of LBGs is not cor- 
rect, as we would be adding the star formation in the cores 
of missed LBGs to the outskirts presumably composed of 
atomic-dominated gas. The true correction would lie some- 
where between the black solid line and the gold long-dashed 
line. We take the molecular gas covering fraction results 
as evidence that the outskirts of LBGs consist of atomic- 
dominated gas, which is consistent with our underlying hy- 
pothesis for this paper. 

5.3. Measurements of the SFR and p„ in the Outskirts of 

LBGs 

In addition to calculating the efficiency of the SFR, we can 
also calculate the mean SFR, (SFR), and p\ in the outskirts 
of LBGs by integrating the rest-frame FUV emission in the 
outer areas. These measurements allow us to calculate the 
metals produced in the outskirts of LBGs in Section 5.4 and 
to put a limit on the tota l p* contributed by DLA gas. Similar 
to I Wolfe & Chenl d2006l) . we assume that DLAs are disk like 
structures, or any other type of gaseous configurations with 
preferred planes of symmetry. We note that while we work 
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Table 2 

SFR, pi, and Metallicity 



01ow a 




R/H b 


SFR 


p* 


p* Corrected 


r 


Z 


[M/H] 


(arcsec) 


(arcsec) 




(Mq yr" 1 ) 


(Mq yr 1 Mpc~ 3 ) 


(Mq yr" 1 Mpc~ 3 ) 




(Z ) 




0.405 


0.765 


10 


0.1 18 ±0.002 


(1.82 ±0.03) x 10~ 3 


(3.75 ±0.03) x 10" 3 


0.084 ±0.001 


0.12 ±0.05 


-0.9 ±0.4 


0.405 


0.765 


too 


0.059 ±0.001 


(0.91 ±0.01) x 10^3 


(1.87 ±0.01) x 10~ 3 


0.084 ±0.001 


0.12±0.05 


-0.9 ±0.4 


0.405 


1.125 


10 


0.1 87 ±0.005 


(2.87 ±0.08) x 10- 3 


(5.91 ±0.08) x Kr 3 


0.126 ±0.002 


0.19 ±0.07 


-0.7 ±0.4 


0.405 


1.125 


too 


0.093 ± 0.002 


(1.44 ±0.04) x 10-3 


(2.95 ±0.04) x 10~ 3 


0.126 ±0.002 


0.19 ±0.07 


-0.7 ±0.4 



Note. — The integrated SFR and p* in the outskirts of LBGs. We integrate from the point where the theoretical models for DLA gas and the LBG data overlap, 
in order to probe the hypothesis that the FUV emission in this region is from in situ star formation in DLA gas associated with the LBGs. 
d Radii from the center of the composite LBG stack. 

^ The thickness of the disk, where R is the radius and H is the scale height. 

c Fraction of p* observed in the outer region of the composite LBG stack divided by the total p* observed. 



this out for DLA gas, the only assumption in our derivation is 
that the outskirts have preferred planes of symmetry such, as 
disks. Even if the outskirts of LBGs are not DLA gas, such an 
assumption is still valid give n the rotation curves m easured 
for high redshift LBGs (e g.. ISchreiberet al 1120091) and the 
predi ctions by simulations (Brooks et al. 2009t lCeverino et ail 
l2OT0h . 

For these calculations, we need S,Jr, the luminosity per unit 
frequency interval per unit area projected perpendicular to the 
plane of the disk. Specifically, we solve for S^r as a function 
of (I°^ s ) averaged over all inclination angles. We find that 

^_ 4^(l + z) 3 (Q 

v ~ MR/H) ' K) 

where R is the radius and H is the scale height of the model 
disks, which holds in the limit R^ H. We calculate S^J: for 
a range in aspect ratios, with R/H values from 10 to 100, 
covering a ra nge from thick disks, as possibly seen at high 
redshift (e.g., Schre iber et al.l l2009). all the way to thin disks 
resembling the Milky Way. We then calculate the mean SFR 
by integrating across the outer region of the LBG stack, 
and find that 



(SFR) = 



8tt(1+z) 3 
C'MR/H) 



2Trd A 2 (C(9))8d9 



(10) 



where dU is the angular diameter distance, 9 is the radius in 
arcseconds, 6>i ow is the minimum radius for the outer region, 
and #hi g h is the maximum radius. I°^(9) comes from the radial 
surface brightness profile from Section 4.3 and depends on 
fi v (9), namely 7° o bs (6») = lO" - 4 ^^ 48 - 6 ). The SFR depends on 

the inclination angles of the disks for a given I°^ s and includes 
a factor of two for averaging over all inclination angles. 

In order to find the (SFR) in the outer regions of LBGs, 
we need to designate a radius at which to start integrating the 
LBG stack, and a comparison of the theoretical model to the 
data in Section 6.2 yields this radius. The (SFR) is indepen- 
dent of the theoretical framework developed later in Section 6 
and does not depend on the efficiency of the gas. It is purely a 
measurement of the star formation occurring in the gas. How- 
ever, it requires a minimum radius to define the beginning of 
the outer region of the LBGs. Specifically, we pick the small- 
est radius that corresponds to the first point in Figure [9] where 
we demonstrate that the smallest radius of atomic-dominated 
gas corresponds to 0."4[3 We note that when we refer to spa- 



tially extended emission star formation throughout the paper, 
we are referring to emission at radii larger than 0."4. We also 
require a second point that we integrate out to, for which we 
use two different values. First, we use the radius of 0."cQ 
which corresponds to the largest radius above 3<r. Second, 
we integrate to 1 ." 1 Q corresponding to the furthest point for 
which we measured py in Figure |4] Table |2]lists the SFRs for 
different combinations of R/H and #high- Since we integrate 
over the radii where the LBG data intersect the theoretical 
models for the DLA gas, the FUV emission from this region 
may be from the in situ star formation in DLA gas associated 
with the LBGs (see Section 7). 

After we have the SFR we can calculate values for the SFR 
per unit comoving volume, p'*, via /9*=SFRxA^lbg/Vudf, and 
they are tabulated in Table [2] However, since p* depends on 
A^lbGi we perform a completeness correction as described in 
Appendix A. 2. The resultant completeness corrected p\ are 
listed in Table [2] While there is a range in the acceptable 
values for both the SFR and p'», we find that the extended 
emission has (SFR) ~O.1M yr" 1 and p> 3 x 10" 3 Mq 
yr" 1 Mpc~ 3 . 

We take this meas ured p* in conju nction with the upper 
limit found in Wolfe & Chen (2006) to calculate the total 
p* fr om neutral atom ic-dominated gas at z ~ 3. Specifi- 
cally, iWolfe & Chenl (120061) constrain p\ for regions in the 
UDF without LBGs, which complements the results from this 
study for regions containing such objects. Together, we con- 
strain all p ossibilities for the star formation from such gas. 
IWolfe & Chenl (120061) place a conservative upper limit on 
p* contributed by DLAs with column densities greater than 
N mia = 2 x 10 20 cirT 2 , finding p\ < 4.0 x 10" 3 Mq yr" 1 Mpc" 3 . 
Combining this with our largest possible value of the com- 
pleteness corrected p'» in Table|2]of p'»=5.91 x 10~ 3 Mq yr -1 
Mpc" 3 , we calculate an upper limit on the total p* contributed 
by DLA gas. We find a conservative upper limit of p'*< 
9.9 x 10" 3 Mq yr" 1 Mpc" 3 , corresponding to ~ 1 0% of the p* 
measu red in the inner regions of LBGs at z ~ 3 (Red dv et alj 
l2008h . 

5.4. Metal Production in the Outskirts of LBGs 

Under the hypothesis that the outskirts of LBGs are com- 
posed of atomic -dominated gas, we can calculate the metals 



10" 3 Mq yr" 1 kpc" 2 (Figure|4). 

13 This corresponds to /xy = 30.5 ±0.5 mag arcsec -2 , or Ssfr : 
0.6) X 10~ 3 Mq yr 1 kpc" 2 (Figure|4). 

14 This corresponds to p,y = 30.6 ± 0.6 mag arcsec -2 , or Sspp> : 



(1.3 ± 



(1.2 ± 



This corresponds to fi v = 28.8 ±0.2 mag arcsec" 2 , or E SFR = (6± 1) X 0.7) X 10 3 Mq yr 1 kpc 2 (Figure[4}. 
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produced due to in situ star formation from z=10toz = 3, 
and compare this to the metals observed in DLAs at z = 3. 
The metal production can be measured from the outskirts of 
the LBG composite since the FUV luminosity is a sensitive 
measure of star formation, since the massive stars produce the 
UV photons as well as the majority of the metals. The co- 
moving density of metals produced is obtained by integrating 
the com oving SFR density (p* ) from the most recent g alaxy 
surveys (Bouwens et al. 2010a b; Reddy & Steidel 2009). We 
note that the resultant metallicities are only valid if the out- 
skirts of LBGs are composed of atomic-dominated gas, as we 
divide by the H I mass density, pm, to obtain the metallicity. 

First, we integrate p\ for all LBGs from z ~ 3 to z ~ 
10 using the p\ values from iBouwens et alj d2010allbh and 
iReddv & Steidell ( 120091) to calculate the total mass of met- 
als pr oduc ed in LBGs by z ~ 3 similar to lPettinil (11999112004 
2006) and lWolfe etafl d2003al) . Specifically, we calculate the 
comoving mass density of stars at z ~ 3 by 



z=io 



dt 



P*,lbg= / p*.LBC-rdz= 1.1 x 1O 8 M Mpc J (11) 
lz=3 dz 



where 



1 



dt 



dz (l + z)H(z) 



(12) 



In order to obtain the comoving mass density of stars in 
the outskirts of LBGs, we multiply this result by the frac- 
tion of p\ observed in the outer region of the composite LBG 
stack compared to the total p\ observed, /, which we list 
in Table [2] We can then calculate the total mass in met- 
als produced by z ~ 3 using the estima ted conversion factor 
Pmetais = (1/64) p'* bv lConti et al.l (120031) . which is a factor of 
1 . 5 lower than the me tal production rate originally estimated 
by lMadau et al.l (119961) . The metallicity of the presumed DLA 
gas is then calculated by dividing by pni at z ~ 3, where we 
use average value of pui over the redshift range 2.4 < z < 3.5 
of (9.0 ±0.1) x 10 7 M Q Mpc" 3 dProchaska & Wolfe 2009h . 
The final metallicities are tabulated in Table [2] in terms of 
the solar metallicity, where Z© = 0.0134 (lAsplund et al.l2009t 
iGrevesse et al.l l2010). The metallicities range from 0.1 2Z© to 
0.19Z Q , similar to DLA metallicities (see Section 7.6). We 
note that / is independent of the disk aspect ratio (R/H), and 
therefore so is the metallicity. 

6. STAR FORMATION RATE EFFICIENCY IN 
NEUTRAL ATOMIC-DOMINATED GAS 

In our search for spatially extended LSB emission around 
LBGs, we aim to further our understanding of the connection 
between the DLA gas studied in absorption and the star for- 
mation need ed to explain the char acteristics of the DLA gas. 
Specifically. IWolfe & Chenl (|2006) searched for isolated LSB 
emission of DLAs, away from known LBGs, in the UDF and 
found conservative upper limits on the SFR per unit comoving 
volume, p\. These limits constrain the in situ SFR efficiency 
of DLAs to be less than 5% of that expected from the KS rela- 
tion. In other words, star formation must occur at much lower 
efficiency in neutral atomic-dominated hydrogen gas at z ~ 3 
than in modern galaxies at z = 0. 

The surface brightness profile of the super-stack of 48 re- 
solved LBGs (Figure reveals the presence of spatially ex- 
tended star formation around LBGs. The latest evidence 



suggests that this star formation is most likely occurring 
in atomic-dominated gas. The most convincing evidence 
is measurements probing the outer disks of local galaxies 
that detect star formation in atomic-dominated h ydrogen gas 
dFumagalli & Gavazzil 12008; Bigie l et alj|2010blab . In addi- 
tion, we find that because the covering fraction of molecular 
gas is insufficient to explain the observed star formation in 
the outskirts of LBGs, the observed emission is likely from 
atomic-dominated gas (see Section 5.2). Throughout the rest 
of this investigation, we work under the hypothesis that the 
observed FUV emission in the outskirts of LBGs is from in 
situ star formation in atomic-dominated gas. In order to quan- 
tify the efficiency of star formation at high redshift in atomic- 
dominated gas, we require a theoretical framework connect- 
ing the observed emission around LBGs to the expectations 
based on known DLA statistics.. 

We develop such a framework in Section 6.1, where we 
combine the column-density distribution function of DLAs 
with the KS relation to construct a model that predicts the 
comoving SFR density per intensity for different SFR effi- 
ciencies of the KS relation. We convert the measured ra- 
dial surface brightness profile from Section 4.3 into this same 
quantity in Section 6.2 and compare it to the model. Through 
this comparison, we obtain an SFR efficiency for each surface 
brightness in the profile which corresponds to both a specific 
radius in the profile and a gas column density via the KS rela- 
tion. As a tool to unde rstand the SFR efficien cies and compare 
our res ults to those of I Wolfe & Chen| ([20 06 ) and Bigi el et al] 
(2010b), and simulations such as Gnedin & Kravt sovl d2010l) . 
we convert our results to fit onto a standard plot of Ssfr ver- 
sus E gas in Section 6.3. The Esfr is determined directly from 
the measured rest-frame FUV flux, and the E gas is determined 
from the S gas , the KS relation, and the SFR efficiency de- 
termined through comparisons of the data with the column- 
density distribution function. 

6. 1 . Theoretical Framework 

We require a theoretical framework for the rest-frame FUV 
emission from DLAs and start with the one developed in 
IWolfe & Chenl d2006l) for the expected emission from DLAs 
in the V band image of the UDF. After summarizing this 
framework, we expand it to explain the observed emission 
around LBGs as a function of radius (and therefore surface 
brightness), taking into account the projection effects of ran- 
domly inclined disks. Our resultant model in §6.1.2 yields 
predictions of the differential p\ per intensity interval ex- 
pected from DLAs for different SFR efficiencies. We then 
compare this model to the data in Section 6.2 to obtain the 
SFR efficiency of the DLA gas. 



6.1.1. Original Framework from 'Wolfe & Chen ( 2006 ) 

The framework developed in IWolfe & Chenl ( 120061) con- 
nects the measured column density distribution function, 
f(Nm,X), the KS relation, and randomly inclined disks to de- 
termine the expected cumulative for DLAs as a function 
of column density and therefore surface brightness. Specif- 
ically, they develop an expression for p\ due to DLAs with 
observed column density greater or equal to N, and we take 
this ex pression directly from Equation (6), in Wolfe & Chen 
(2006), namely 



dN'j(N')X SFR (N') . (13) 
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Here Hq is the Hubble constant, c is the speed of light, N max 
is the maximum observed column density for DLAs (10 22 
cm" 2 ), and J(N') is 



J(N') = 



min{No,N ) 



dN ± g(N ± ,X) 



N± 



N' 



(14) 



Here X(z) is the absorption distance and g(N±,X) is the in- 
trinsic column-density distribution of the disk for which the 
maximum value of N±, the H I column density perpendicu- 
lar to the disk, is A^o- g(N±,X) is related to the observed H I 
column-density distribution function fm(N,X) by 



fm(N,X) = 



min(No,N) 



dN±g(N±,X)(Nl/N*) (15) 



( Fall &Peilll993t IWolfe et alJl995T) . 

A potential problem with using f(Nm,X) in the expression 
for pi is that the measurements of f{Nm,X) originate from 
absorption-line meas urements that sample scales of ~ 1 pc 
dLanzetta et all [2002). On the other hand, the KS relation 
is est ablished on scales exceeding 0.3 kpc (Ke nnicutt et al.l 
120071) . This is not an issue, however, because f(Nm,X) typ- 
ically depends on over 50 measurements per column-density 
bin and is therefore a statistical average over probed areas that 
exceed a few kpc 2 (see IWolfe &"Chenll2006l) . 

6.1.2. New Differential Approach for LBG Outskirts 

The framework developed in Wol fe & Chenl (120061) was ap- 
propriate for connecting the upper limit measurements of pi 
from DLAs above a limiting column density and therefore 
surface brightness, to model predictions based on f{Nni,X) 
and the KS relation. It does not, however, work in the 
present context of positive detections over a range of surface 
brightnesses. For this we require a differential expression for 
pi, rat her than a cumulative version used by Wolfe & Chenl 
(2006). Specifically, assuming that LBGs are at the cen- 
ter of DLAs, we wish to predict the rest-frame FUV emis- 
sion of DLAs for a range of efficiencies of star formation 
in such a way that we can distinguish between possible dif- 
ferent efficiencies for each surface brightness interval. We 
find that dpl/dN accomplishes this by yielding unique non- 
overlapping predictions for each efficiency. Each differential 
interval of pi represents a ring around the LBGs correspond- 
ing to a surface brightness and a solid angle interval subtended 
by each ring. This surface brightness corresponds to the col- 
umn density of gas corresponding to some radius in the radial 
surface brightness profile and is responsible for the emission 
covering that area on the sky. If this gas is neutral atomic- 
dominated H I gas, then we can predict the expected dpl/dN 
using f(Nni,X) and the KS relation. 

Specifically, to obtain dpl/dN, we differentiate Equation 
(\3[ with respect to N. To do so, we need g(N±,X), since 
Equation ( fT3l depends on Equation (fT4l . which depends on 
g(N± , X) . We find g(N± , X ) from f(Nm ,X)to obtain a general 
form of the equivalent double power-law fit for g(N± , X) when 
N± < No using Equation (fTBT l. We find that 



(AT \ a 
at) ' n±<n ° 



(16) 



g(N±,X)=0 forN± > Nq. We now differentiate Equation ( fT3l 
with respect to N to get 



dpi hfm fHo\fKk 3 



a + 3 



(3 + 2 + a J XN 2 



where 



h(N) : 



N P + 2 +a _ N P+2 +a . N<N ^ 



N 



/3+2+a 







■NT a ;N>N 



(17) 



(18) 



In the case of a =a 3 =- 2.00±0.05 for N < Nq 

dProchaska & Wolfd 120091) . this reduces to 

g(Nj_,X)=k3(N±/NoT 2 at N± < No and g(N±,X)=0 for 
N± > No- Also, Equations [T7land[T8lreduce to, 



d -£l= K N)(^ 
dN V c 



-2) (19) 



where 



h(N) = 



N P -N p mm ;N<N , 
N?-K; n ;N>N , 



(a = -2) 



(20) 



We note that the expression for dpl/dN is independent of 014. 

We would like to compare dpl/dN to the observations, 
however, we cannot measure dpl/dN directly. On the other 
hand, we can measure d pi / d (I°^ s ) , which is easily derived 
from dp* /dN. Specifically, we find 



dpi 

d(m 



dpi 
~dN 



dT, 



SFR 



dN 



( dT, 



SFR 



Since 



and 



dE 



SFR 



dN 



K/3N 13 - 



afEsFR 4^(l+z) 3 /3 



d(I%* 
we therefore find that 



C 



(21) 



(22) 



(23) 



a + 3 
13 + 2 + a 



( No a 
\N' M 
(24) 



d{^) """\c)\ 
which for a = -2 reduces to 

dW) im \~)\~^^)\^)' ( ) 

where h(N) is the same as in Equations ( fT8l and d20b . Since Af 
is related to (I°^ s ) through Equation © and the KS relation, 
we find 



N: 



b ■ ( 4n(\+z?Wc 

robs 



(26) 



where £3 and a are the same as in Equation ||7), and thus surface brightness. 



and therefore d pi / d is a unique function of (/* s ), and 
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Figure 7. Solid blue curve is the surface brightness (fiy ) versus the dif- 
ferential comoving SFR density per intensity Ap'*/A{/°[j s ) predicted for 
Kennicutt-Schmidt relation with K = K^ am . The red curves depict what 
Ap'*/A(/°jj s ) looks like for different values of the normalization constant 

K, where ^Kenn=(2.5±0.5)x 1(T 4 Mq yr 1 kpc~ 2 <Kennicutl|[T9983bT) . The 
curves in this figure predict the amount of star formation that should be ob- 
served around LBG s for different efficiencies of star formation from Equa- 
tions <25l and (26). The dotted line is for K = 0.1 X K K enn , the dashed line 
is for K = 0.05 X Kkenii! the the dot-dashed line is for K = 0.02 X XKenn, and 
the tripple-dot-dashed line is for K = 0.01 X Xxenn- The range of surface 
brightnesses corresponds to 5 X 10 20 cm -2 < N < 1 X 10 22 cirT 2 . 



The resulting predictions for the surface brightness, py, 
versus the differential comoving SFR density per intensity, 
d p'* I d (I°^ s ) , is depicted by the blue curve in Figure|7] The red 
curves in this figure depict dp\/d{I^) for different values of 
the normalizatio n constant K, whe re ^Kenn=(2.5±0.5)x 10~ 4 
M Q yr" 1 kpc" 2 dKennicutu Il998allbl) . First, we note that 
dp*/d(I^ s ) for a given K decreases as py decreases because 
of the decreasing population of high column-density DLAs 
(i.e. DLAs with higher surface brightnesses). Second, we 
note that both p* and are linearly proportional to K, and 
therefore K cancels out in the x-direction, leaving only the ob- 
served pv to vary with K in the y-direction. We plot pv on the 
>'-axis to conceptually facilitate the conversion of these results 
later in the paper and note that d p\ / d (1°^) is a function of N 
and therefore py. The curves in this figure predict the amount 
of star formation that should be observed around LBGs for 
different efficiencies of star formation, which will be com- 
pared to the data in Section 6.2. We plot d p* / d (/? s ) for 
the range in py that corresponds to N m i n < N < N m!lx , where 
N mia = 5 x 10 20 cirr 2 , and AUx = 1 x 10 22 cm -2 . The value 
of A^ m in is lower than the range of threshold column densi- 
ties of N" n observed for nearby galaxies (lKennicutt| [l998b). 
but at a column density high enough such that we may start 
to see star formation occur. Also, recent results probing 
in the outer disks of nearby galaxies observe star formation 
at low column densities in atomic-dominated h ydrogen gas 
dFumagalli & Gavazzil 12001 iBigiel et al.l l2010bllah . Regard- 
less, our measurements do not probe column densities down 
to this level, so the exact value for N m i n does not affect the 
results. 



obs\ 



6.2. Stacking Randomly Inclined Disks 

We wish to compare these theoretical values of dp\jd(J\ 
predicted for DLA gas forming stars according to the KS re- 
lation, to empirical measurements of Ssfr for our LBG sam- 



ple. To do this, we require a method to convert the radial sur- 
face brightness profile in Figure|4]into py versus d p* / d . 
For a given ring corresponding to a point in the radial surface 
brightness profile and covering an area AA we can calculate 
Ap* from the measured intensity (7J? S ). Ap'» is similar to 
the differential p* mentioned above, and is calculated from 
the flux measured in an annular ring at some radius from the 
center of the LBGs. Specifically, 



A/9'* : 



AL^Vlbg 

C'Vudf 



(27) 



where C' is the conversion factor from FUV radiation to 
SFRB C = 8 x 10 27 erg s" 1 Hz" 1 (M yr" 1 )" 1 , AL U is the lu- 
minosity per unit frequency interval for a ring with area AA, 
A^lbg is the number of z ~ 3 LBGs in the UDF, and Vudf is 
the comoving volume of the UDF. As discussed in Section 
3.1, we use a comoving volume of 26436 Mpc 3 . We recog- 
nize that A^lbg is dependent on the depth of images available 
to make selections, and discuss this completeness issue in Ap- 
pendix A. 3. The value of Ap'„ depends on AL„, and therefore 
on the inclination angle i for a given measured intensity, in the 
case of planes of preferred symmetry. 

In order to determine Ap'*, we average over all inclination 
angles in our determination of AL^, which depends on 
described in Section 5.2. Specifically, we use Equation (O in 
conjunction with AL V = AAj^E^ to find Ap'*, where AA^ 
is the area parallel to the plane of the disk. First, we rewrite 
AA^ in terms of AA, the projection of AAj^ perpendicular 
to the line of sight, and find that averaging over all inclination 
angles yields AA = 2AA. We can then rewrite AA^ in terms 
of Af2, the solid angle subtended by one of the rings from the 
surface brightness profile, and AA, yielding AA 
Using this relation, we find 



2AfW 2 . 



Ap; 



87r(l+z)V A ^ L B G AF° bs (r) 
\n(R/H)C'V VDF 



(28) 



where AF° bs (r) is the observed integrated flux in a ring as a 
function of the radius r. We then calculate the change in in- 
tensity from one ring to the next, A(/* s ), to get Ap;/A(/* s ) 
by measuring the intensity change across each ring by taking 
the difference between values of the intensity on either side 
of each point and dividing by two. In the case that A(/°^ s ) as 
calculated above is negative, we take the change in intensity 
over a larger interval. 

Figure[8]shows a comparison of the theoretical model from 
Section 5.2.2 and the measured values from the radial sur- 
face brightness profile. The measurements are for a range in 
aspect ratios, with R/H values from 10 to 100, similar to Sec- 
tion 5.2. We display the data in two complementary ways. 
First, the green diamonds depict results assuming the average 
value of the possible aspect ratios, with the error bars reflect- 
ing the measurement uncertainties (including the uncertainty 
due to the variance in the image composite) in order to portray 
the precision of our measurements. Second, we show a filled 
region, where the gray represent results for the full range in 
aspect ratios and the gold represents the la uncertainties on 
top of that range. This shows the region that is acceptable 
for each of those points. In both portrayals, we only include 
uncertainties of the aspect ratios, the variance due to stacking 

15 C' is in different units than the same factor C in Equation 
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Figure 8. Surface brightness versus the differential comoving SFR den- 
sity per intensity (Ap'»/A{/°[j s )), comparing the measured emission in the 
outskirts of LBGs to the predicted levels for different SFR efficiencies. The 
blue and red lines are the predictions from Figure [jj The filled gray region 
represents observed emission in the outskirts of LBGs for a range in aspect 
ratios, with R/H values ranging from 10 to 100, and the filled gold region is 
its lcr uncertainty. The green diamonds are the average value of the possi- 
ble aspect ratios, with the error bars reflecting the measurement uncertainties 
(including the uncertainty due to the variance in the image composite). 

different LBGs, and the measurement uncertainties, and do 
not include uncertainties in the FUV light to SFR conversion 
factor (C) from Equation (0]), or any other such systematic un- 
certainties. We truncate the data at a radius of ~ 0.8 arcsec 
(~ 6 kpc) corresponding to a 3<r cut to include only measure- 
ments with high signal to noise. We note that in calculating 
the S/N values, we include the uncertainties due to the vari- 
ance of objects in the composite stack. The data beyond a ra- 
dius of ~ 0.8 arcsec are plotted in gray, and while they yield 
similar results, they are not included due to their low S/N. 

The results shown in Figure [8]do not yet include complete- 
ness corrections, which we describe in Appendix A. 3. We 
present the completeness corrected comparison between the 
theoretical models of d p* / d (I°^ s ) for DLAs with measured 
values of Ap* / A(I°^ S ) in Figure|9] This shows that, under the 
hypothesis that the observed extended FUV emission comes 
from in situ star formation of DLA gas, the DLAs have an 
SFR efficiency at z ~ 3 significantly lower than that of local 
galaxies[3- In fact, the Kennicutt parameter K needs to be re- 
duced by a factor of 10-50 below the local value, K = K^ mn . 
The values of Ap'*/ A(I°^ S ) that intersect the predictions of 
the theoretical models for N = 5 x 10 20 — 1 x 10 22 are black 
crosses and have S/N values ranging from ~17 to ^3 sug- 
gesting that the measurements are robust. These points cor- 
respond to radii of ~ 0.4-0.8 arcsec (~ 3-6 kpc). The point 
with the largest value of Ap*/A(I°^ s ) of — 17.4 in Figure|9] 
seems to deviate from what otherwise would be a clear trend. 
This point corresponds to the point at a radius of 0.72 arcsec 
in Figure |4] which also differs slightly from the general de- 
creasing trend in py However, it is consistent within their 
uncertainties for the surface brightness profile, and we are not 
concerned about it. All the points at radii larger than ~ 0.8 
arcsec also intersect the theoretical models with similar ef- 
ficiencies, but are not included as they have lower S/N. The 
values of K vary for each data point and are not constant for a 

16 We remind the reader that if the working hypothesis is not correct, then 
the results of Wolfe & Chen 1 2006) already constrain the SFR efficiency. 
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Figure 9. Same as Figure[8] but corrected for completeness. The points that 
overlap with the theoretical models are now black crosses, to emphasize the 
points that will be used for the rest of this work. These points correspond 
to radii of 0."405 through 0."765. We omit the low S/N (< 3a) gray data 
points from Figure [8] By comparing the measurements in this plot to the 
predictions, we get the SFR efficiency for each surface brightness. 

given pv- These tantalizing results will be discussed in Sec- 
tion 7, and we consider the effects of stacking different sam- 
ples of LBGs in Appendix B. 

6.2. 1 . Variations in the KS Relation Slope f} 

The SFR efficiencies can also be decreased by lowering the 
slope (3 of the KS relation, while keeping K = .^Kenn- The value 
of P in the l iterature ranges from /} ~ 1 .0 dBigiel et alj|2008l) 
to /3 ~ 1.7 (Bo uche et al.l 12007b . Increasing the value of j3 
would increase the SFR efficiency, so we do not consider that 
here. On the other hand, lower values of j3 would decrease the 
SFR efficiency and there ar e physical motivations to conside r 
/3 values as low as 1.0 (e.g., Elmegreen 2002; Kravtsov 2003). 
However, even reducing j3 to values as low as 0.6 does not 
reduce the SFR efficiency enough to match our observations, 
and there are no physical motivations nor data to justify values 
of j3 lower than 0.6. Lastly, decreasing the value of j3 not only 
decreases the SFR efficiency, but it also decreases the value of 
dp'*/d(I°^ s ) (see Equation d25ll). While decreasing j3 would 
move the blue curve in Figure [9] down, it also moves it to the 
right and does not help the models match our observations. 
Therefore, we focus on other mechanisms for reducing the 
SFR efficiencies by varying the value of the parameter K. 

6.3. The Kennicutt Schmidt Relation for Atomic-dominated 
Gas at High Redshift 

As a tool to understand the low SFR e fficiencies, and in or- 
der to compare our resu lts to those of Wolf e & Chenl (120061) 
and Biaiel et al. (2010b), and simulations such as those by 
Gnedi n & Krav tsov (2010), we translate the result from Sec- 
tion 6.2 and Figure [9] to a common set of parameters to obtain 
a plot of Ssfr versus E gas . We derive this conversion for the 
emission in the outskirts of LBGs in Sect i on 6.3 .1, and for 
the DLA upper limits from I Wolfe & Chenl (l2006h in Section 
6.3.2. 

6.3. 1 . Converting the LBG Results to Esfr versus E gas 

We calculate Esfr directly using Equation (0|, where the 
average intensity is obtained from the radial surface bright- 
ness profile of the composite LBG stack in rings as explained 
in Section 4.3 (Figure|4|. As p v algebraically increases with 
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Figure 10. Star formation rate per unit area (Xsfr) versus gas density (S gas ). The dashed line represents the KS relation with K = XKcnn=(2.5±0.5)x 10 
Mq yr _1 kpc~ 2 , while the triple dot-dashed line is for K = 0.1 X ^Kenn- The gray filled region, the gold filled region, and the black cros ses represent the same 
data as in Figure H The green data points corresponding to upper limits derived for DLAs without central bulges of star formation from Wolfe & Chen 1 2006) 
converted to work with this plot. Since DLA sizes are not known, the upper limits are derived for angular diameters <9kern=4", 2", 1", 0."6, 0."5, 0."4, 0."3, 
and 0."25, from left to right. The data all fall at or below 10% of the KS relation, showing a lower SFR efficiency than predicted. 



increasing radius, Ssfr will generally decrease with increas- 
ing radius. We compute S gas by first inferring the value of K 
for each data point that intersects the theoretical curves in the 
Hv versus Ap'*/A(/°^ s ) plane, since each theoretical curve is 
parameterized by a fixed value of K. We precisely find the 
corresponding efficiency by calculating a grid of models with 
K varying by 0.001. We then use this value of K and the KS 
relation (Equation (fJI) to calculate E gas for each value of Ssfr 
inferred from the measured fly . 

In short, we directly measure Ssfr through the emitted 
FUV radiation, and then calculate S gas for the correspond- 
ing K that matches the DLA model. We note that Ssfr is a 
direct measurement, while S gas comes from the DLA model 
which is based on the column-density distribution function 
of DLA gas. The result is shown in Figure [TUl which plots 
Ssfr versus S gas . We truncate the plot to include only data 
that overlaps with observed DLA gas densities, namely N < 
1 x 10 22 cm" 2 , where N is calculated from Equation (01 using 
the K value determined for that /j,y. The dashed line black 
represents the KS relation with K = /T K enn=(2.5±0.5)x 10" 4 
Mq yr _1 kpc~ 2 , while the pink triple dot-dashed line repre- 
sents K = 0.1 x ^Kenn- The gray filled area with the la un- 
certainty and the black points represent the same data as in 
Figure|9] The green upper limits will be described in Section 
6.3.2. 

The LBG outskirts in Figure[l0]clearly have lower SFR ef- 
ficiencies than predicted by the KS relation. In addition, they 
appear to follow a power law that is steeper than the KS rela- 
tion at low redshift. However, there is a large scatter caused by 
the uncertainty introduced by the sky-subtraction uncertainty 
(see Section 4.2), and the sample variance due to stacking dif- 
ferent objects (see Section 4.3). We are therefore cautious 
about fitting a power law to this data by itself, and investigate 
this trend further in Section 7.3.1. 

6.3.2. Converting DLA Data to Esfr versus E gas 

To convert the DLA points from I Wolfe & Chenl ( 120061) . we 
use the same idea as that of the LBG data, except that in this 
case we do not have detected star formation, and therefore 



we use upper limits. Also, rather than using the new frame- 
work developed above, we use the formalis m developed for 
DLAs without central bulges of star formation (IWolfe & Chenl 
Hj(p. We start with Figure 7 o flWolfe & Chenl (120061) . which 
basically plots the SFR density due to star formation in neu- 
tral atomic-dominated hydrogen gas (i.e., DLAs) with column 
densities greater than N (p'*(> AO) versus M v Pl Similar to our 
FigureEl Figure 7 of I Wolfe & Chenl d2006l) has DLA models 
with different Kennicutt parameters K. We follow the same 
technique of using the intersection of these models with the 
data to find K values for each /iy. However, in this case the 
values of \M/ do not correspond to detected surface bright- 
nesses of LBGs, but rather correspond to threshold surface 
brightnesses, i.e., the lowest values of (I°^ s ) that would be 
measured for a DLA of a given angular diameter. Therefore, 
we cannot calculate Ssfr the same way as in Section 6.3.1. 

To determine Ssfr, we first calculate an effective minimum 
column density, AW, corresponding to the threshold surface 
brightness of IWolfe & Chenl d2006l) . We do this using the 
K values from the intersection points, the thresholds /iy, and 
Equation[3] We then calculate Ssfr, the average over all pos- 
sible column densities above AW, as 



Ssfr = (S SF r(> AW)) = 



C / 0'»(>A e ff) 



N -J(N')dN' 



(29) 



where J{N) is the integral in Equation ( TBI ). We then calculate 
S gas using these Ssfr, the above K values, and Equation ©. 
The resulting values are overlaid on the LBG results in Figure 
[TOl We emphasize that the two results are independent tests. 
The DLA upper limits put constraints on the KS relation in 
the case of no central bulge of star formation, while the LBG 
outskirts data put constraints on in situ star formation in DLAs 
associated with LBGs. Together, these results show that the 
SFR efficiency in diffuse atomic-dominated gas at z ~ 3 is less 
efficient than predicted by the KS relation for local galaxies. 

17 The .v-axis of the plot is actually Esfr, but Ssfr corresponds to a value 
of nv, which is easier to understand and makes more sense in this context. 
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Figure 11. S tar formation rate p er unit area (Ssfr) versus gas density (S gas ) for both the data as shown in Figure[lO] and the KS relation for simulated galaxies 
at z = 3 from I Gnedin et al. 2009). The left panel includes gas of all metallicities while the right panel considers only the gas with metallicities below 0. IZq . The 
solid red and blue lines are the mean relation for the total neutral-gas surface density (atomic and molecular) and the hatched area shows the rms scatter around 
the mean. The dotted purple line shows the mean KS relation only consid ering the ato mic gas, while the dot-dashed cyan line shows it for the molecular gas. The 
long-dashed line is the Kennicutt relation for local (z ~ 0) galaxies I Kennicutt 1998a), while the triple-dot-dashed line is for K = 0.1 X ^Kenn- The gray and gold 
lines and the black and green points are the data and are the same as in Figure [Tol 



7. DISCUSSION OF RESULTS 

We present, for the first time, evidence for LSB emission 
around LBGs on spatial scales large compared to the LBG 
cores due to in situ star formation in gas associated with LBGs 
(see Figure|4|. Using the theoretical framework developed in 
Section 6.1.2, we model this emission as the average in situ 
star formation in atomic-dominated gas surrounding LBGs at 
Z ~ 3. Since most of the atomic-dominated gas at high red- 
shift is in DLAs, we assume that this is DLA gas. We find 
that under this hypothesis, the efficiency of star formation in 
the atomic-dominated gas is significantly lower than what is 
expected for predictions by the KS relation. This is clearly 
evident in Figure [10] which compares the local KS relation 
directly to the observations. The efficiency of in situ star for- 
mation in atomic-dominated gas appears to be a factor of 10- 
50 lower than that of local galaxies that follow KS relation. In 
addition, the SFR efficiency of atomic-dominated gas around 
galaxies wit hout bright core s is con strained by the thresholds 
measur ed by Wolfe &~Chenl (|2006), who impose similar effi- 
ciencies on the in situ star formation in DLAs, as shown by the 
green upper limits in Figure [10] Therefore, together with the 
results from Wolfe & Chenl (120061) . we constrain the SFR effi- 
ciency of all neutral atomic-dominated hydrogen gas (DLAs) 
atz~3. 

There are multiple possible effects contributing to the ob- 
served lower SFR efficiencies, including a higher background 
radiation field at high redshift, low-metallicity of DLAs, and 
the role of molecular versus atomic hydrogen in star for- 
mation. We start the discussion by comparing our results 
to the models of low-me tallicity high-redshift galaxies by 
iGnedin & Kravtso"vl (120101) in Section 7.1. We then consider 
the roles of molecular and atomic-dominated gas in the KS 
relation in the context of the saturation of atomic dominated 
in Section 7.2. We investigate whether there is a variation of 
the KS relation with redshift in Section 7.3, and compare our 
results at high redshift with those in the outer disks of local 
galaxies in Section 7.3.1. We then consider the effect that a 
bimodality of the DLA population would have on the results 
in Section 7.4, and a caveat on the results based on DLA sizes 
in Section 7.5. Lastly, we address the "Missing Metals" prob- 
lem of DLAs in Section 7.6. 



7. 1 . Models of the Kennicutt Schmidt Relation in 
High-redshift Galaxies 

Here we discuss two effects contributing to the observed 
lower efficienci es: First, at higher reds hifts the background 
radiation field ( Haardt & Madau 1996) is stronger, yield- 
ing a higher UV-flux environment. This photodissociates 
the molecular hydrogen (H2) content of the gas, raising the 
threshold for the gas to become molecular, therefore requir- 
ing higher gas densities to form stars. Second, the metallic - 
ity of DLAs at high redshift is con s iderably lower than so- 
lar dPettini et all 11994 [19951 11997L l2002at iProchaska etafl 
120031k and therefore has a lower dust content, which is needed 
to form molecular hydrogen and to shield the gas from pho- 
todissociating radiation. 

Recent t heor etical wo rk (iKrumholz et alJ 12008 ; 
Gnedin etlfl____|200l IKrumholz et al.1 l2009blal : 
Gnedin & Krav tsovl 120101 12011b suggests that the most 
important part in determining the amplitude and slope of 
the KS rela tion is the dust abundance , and therefore the 
metallicity. Gnedi n & Kravtsovl (2010) investigate the KS 
relation at high redshift usin g their metallicity dependent 
model of molecular hydrogen (iGnedin et al.l l2009). They find 
that while the higher UV flux does lower the SFR, it also 
lowers the surface density of the neutral gas, leaving the KS 
relation mostly unaffected. However, they find that the lower 
metallicity, and therefore lower dust-to-gas ratio, causes a 
steepening and lower amplitude in the KS relation. This 
yields lower molecular gas fractions, which in turn reduces 
SFR for a given gas surface density, S gas . This would yield 
lower observed SFR efficiencies similar to what is measured 
in this study. 

We compare our results with the z = 3 model KS relation 
(Gnedin & Kravtsov 2010) in Figure QT] with a plot showing 
our data overlaid with the model results, and find that their 
predictions are c onsistent with our findings. This plot is sim- 
ilar to Figure 3 in lGnedin & Kravtsovl d20 10). as we provided 
them with our preliminary results for comparison. However, 
there was previously an error in our implementation of the 
theoretical framework, yielding slightly different results than 
presented here. Also, here we split the figure into two panels, 
with the left panel including gas of all metallicities for galax- 
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ies at z = 3, and the right panel including only the gas with 
metallicities below 0.1Z Q . Each plot includes the same points 
as shown in Figure[lO] In addition, the plots show the KS rela- 
tion for the simulated galaxies at z = 3 for the total neutral gas 
(atomic and molecular), with the solid lines showing the mean 
value and the hatched area showing the rms scatter around the 
mean. The dotted purple line shows the mean KS relation only 
considering the atomic gas, while the dot-dashed cyan line 
shows it for the molecular gas. The dashed black line is the 
best-fit relation of iKennicuttl (1 1 998al) for z ~ galaxies, while 
the triple dot-dashed line is for K = 0.1 x ^Kenn- The blue line 
in the right panel is closer to the ran ge of observed metallic i- 
ties observed in DLAs of — 0.04Z Q dProchaska et al.ll2003l ; 
M. Rafelski et al. 2011 in preparation). The use of a metal- 
licity cut of 0.1Z Q is reasonable, since the mass-weighted and 
volume-weighted metallicity of the atomic gas is ~ 0.02Z Q 
and ~ 0.03Z Q respectively (N. Gnedin 2010, private com- 
munication), matching the observed DLA metallicities nicely. 
In addition, although the dispersion of DLA metallicities is 
large, the majority of DLAs have metallicities below 0.1 Z©, 
making this a good choice for a cut to compare to DLA gas. 

The spatially extended emission around LBGs is consistent 
with both the total and low-metallicity gas models in Figure 
[TT1 While the uncertainties in both the model and the data are 
large, taking the results at face value, we can gain insights into 
the nature of the gas reservoirs around LBGs. The data are a 
better match to the model with the metallicity cut of 0.1Z Q , 
coinciding with the mean relation for this model. On the other 
hand, the emission from gas around LBGs is also consistent 
within the la uncertainties of the model including gas with 
all metallicities for the part of the rms scatter below the mean. 
This lower part of the hatched area in the left panel of Figure 
[TTIrepresents th e gas at the lower en d of the metallicity distri- 
bution shown in lGnedin et al.1 (120091) . having a mass-weighted 
and volume-weighted metallicity of the atomic gas of ~ 0.26 
Z© (N. Gnedin 2010, private communication), which is con- 
sistent with that o f z = 3 galaxies having an average metallic - 
ity of - 0.25 Z (IShaplev et al.ll200l iMamiucci et al.ll2009h . 
Since the observations fall below the mean model for the 
higher metallicity model, and coincide with the model for the 
lower metallicity model, we conclude that the metallicity of 
the gas is most likely around the mean metallicity of the low- 
metallicity model, ~ O.O4Z , and is definitely below ~ 0.26 
Zq, the metallicity of the model including gas of all metallici- 
ties The results therefore imply that the gas in the outskirts of 
LBGs has lower metallicities than the LBG cores. In fact, the 
observed SFRs and the implied metallicities from the models 
are fully consistent with the outer regions of LBGs consisting 
of DLA gas. 

7.2. The Roles of Molecular and Atomic-dominated Gas in 
the Kennicutt Schmidt Relation 

The reduction in SFR efficiency will also be affected by 
the competing roles played by atomic and molecular gas in 
the KS relation. There has been some debate as to whether 
the KS relation should include both atomic and molecular 
gas. All stars are believed to form from molecular gas, and 
some workers argue that Ssfr correlates better with molecu- 
lar gas than atomic gas dWong & Blita l2002: Ke nnicutt et al.l 
120071: iBigiel et al.ll2008l) . On the other hand, other observa- 
tions show a cle ar correlation of the total gas surface density , 
S gas , with Ssfb (lKennicutj|1989L [T998al: ISchuster et a"T] |2007t 
Crosthwaite & Turner 20 Q3), and the ato mic gas surface den- 
sity. Shi, with Ssfb (IBigiel et al.ll2010bl) . 



These differences can be understood in the context of the 
saturation of atomic-dominated gas at high column densities. 
Above a specified threshold of Shi, the atomic gas is con- 
verted into molecular gas and no longer correlates with Ssfr. 
This saturation of atomic hydrogen gas above a th reshold sur- 
face d ensity is clearly observed (see Figure 8 of Bigi el et al.1 
2008) for local galaxies at z = 0, and occurs at surface densi- 
ties of ~ 1OM pc -2 (IWong & Bhtzll2002t IBigiel et al.ll2008l) . 
However, below this threshold surface density, there is a cor- 
relation of Shi and Esfr- This is most clearly se en in the re- 
sults an alyzing the outer disks of nearby galaxies (Big iel et al.l 
2010b), where a clear correlation of Shi with Ssfr is ob- 
served. In fact, they find that the key regulating quantity 
for star formation in outer disks is the column density of 
atomic gas. Further evidence is seen in the outskirts of M83, 
where the d istribution of FU V flux again follows the Ehi 
(Bigie l etal.ll2010al) . In fact, IBigiel et all (l2010al) find that 
in the outskirts of M83, massive star formation proceeds al- 
most everywhere H I is observed. While these outer disks 
must contain some molecular gas in order for star formation 
to occur, they are nonetheless dominated by atomic gas with 
surface den s ities l ower than the saturation threshold seen in 
Bigiel etal.1 (120081) . 

The saturation of atomic gas above a t hreshold surface 
density is investigated in theoretical mo dels (iKrum holz et al. 
l2009at lGnedin & KravtsovllMollMTl) . which reproduce the 
saturation threshold for atomic gas in local galaxies. More- 
over, the authors find that the threshold surface density for 
saturation varies with metallicity, where lowe r metallicity sys- 
tems h ave higher thresholds (see Figure 4 of Krumh olz et al.1 
2009a ). In addition, the simulations of z = 3 galaxies by 
Gnedi n & Kravtsovl (120101) show an increased saturation sur- 
face density of atomic gas of ~ 50M Q pc" 2 , as seen in 
the purple dotted line in Figure QT| The line clearly satu- 
rates, with no clear relation between Shi and Ssfr above 
~ 50M© pc" 2 , but with a clear relation below this density. 
Since z = 3 galaxies are typically LBGs with metallici ties of 
~ 0.25Z Q (Sha plev et all 120031: iMannucci etail 120091). thi s 
is consistent with predictions from IKru mholz et al. (2009a). 
The model with only lower metallicity gas (right panel of 
Figure [TTJ saturates at even h igher atomic gas surface den- 
sities (Gnedin & Kravtsovl 120101) . continuing the relation of 
the threshold saturation of Shi with metallicity. It appears 
that Shi tracks S gas to large values of Ssfr, and no saturation 
occurs in S H i until ~ 200-300 M pc~ 2 . 

Comparison of the saturation thresholds found by 
iGnedin & Rravtsov (2010) with our data in the right panel 
of Figure [TT]reveals that the data are below this threshold for 
saturation of atomic gas for all values of Shi for the model 
with the O.IZq cut, which most closely matches our results. 
We therefore conclude that the measured decrease in SFR ef- 
ficiency is not due to the saturation of atomic-dominated gas. 
In fact, our data confirm that the atomic-dominated gas does 
not saturate for Shi of at least > 100M© pc" 2 . This number is 
larger than predicted for LBG metallicities, yet smaller than 
predicted for DLA metallicities. This suggests that the aver- 
age metallicity of the gas in the outskirts of LBGs is between 
0.1Z© and 0.25Z Q . This is consistent with the metallicity es- 
timate made in Section 5.4 of 0.12Z Q to 0.19Z Q based on the 
metal production rate of star formation. 

While molecular gas is very likely needed to form stars, 
it is als o clear from results for local outer disks (Bigi el et al.l 
2010b} and our results at high redshift that star formation can 
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be observed even when molecular gas does not dominate az- 
imuthal averages in rings of galactic dimensions. We do not 
suggest that stars are forming from purely atomic gas, but 
rather that in the presence of high density atomic gas, there 
are sufficient molecules to cool the gas such that star forma- 
tion can be initiated by gravitational collapse. We predict that 
more sensitive measurements of molecular gas in local outer 
disks would result in the dete ction of molecu l ar gas . On the 
other hand, recent models by lLow & Gloved (120101) suggest 
that molecular hydrogen may not be the cause of star forma- 
tion, but rather a consequence of the star formation. In this 
scenario, CO traces dense gas that is already gravitationally 
unstable and would probably form stars regardless. More in- 
vestigations are needed to better understand what is funda- 
mentally required for star formation, and our work provides 
some observational constraints for such work. 

We note that even though large amounts of m o lecula r 
gas are not observed in DLAs (iCurran et al.1 120031 12004) . 
the molecular gas has a very small covering fraction, and 
therefore is unlikely to be seen along random sight lines 
(Zwaan & Prochaska 12006b . We also note that Bigi el et all 
d20 1 Obi) find that the FUV emission reflects the recently 
formed stars without large biases from external extinction. 
We similarly do not expect much extinction in the out er parts 
of the LBGs, as DLAs have low dust-to-gas ratios (Pettini 
2004; Frank & Peroux 2010). Therefore, unless the gas in 
the outskirts of LBGs is due to star formation in molecular- 
dominated gas, we do not expect large extinction corrections 
to be necessary. 

7.3. Is There a Variation of the Kennicutt Schmidt Relation 
with Redshift? 

Recent studies of star-forming galaxies at high redshift have 
sugges ted that the KS relation does not vary with redshift 
(Bouc he et all 120071: iTacconi et ai]l2010t iDaddi et al J 120101: 
Genz elet al.1 120 lOL These studies do a careful job of com- 
paring Ssfr and the molecular gas surface density, £h 2 > of 
high- and low-redshift systems, and find that the galaxies fit 
a single KS relation. This is starkly different than our finding 
a lower SFR efficiency at high redshift. Furthermore, these 
studies differ from simulations o f high-redshift galaxies tha t 
do find a reduced SFR efficiency (Gnedin & Rravtsov 2010). 

These differences may be due to comparisons of different 
types of gas. In our results, we consider atomic-dominated 
gas as found in DLAs, while the other observational studies 
are focused on molecular-dominated gas with galaxies hav- 
ing high molecular fractions a nd higher metal abundances. 
Similarly, Bigi el et al.1 (HpiOb) also found a lower SFR ef- 
ficiency in the outskirts of otherwise normal local galaxies 
when probing atomic-dominated gas. There are two differ- 
ent possible scenarios that explain the results. First, the KS 
relation for atomic-dominated gas follows a different KS re- 
lation than the KS relation for molecular-dominated gas. This 
possibility would explain the observed lower efficiencies of 
star formation in (1) the outskirts of z ~ 3 LBGs as measured 
he re, (2) the DLAs with out star-forming bulges as measured 
bv I Wolfe & Chenl (120061) and (3) the o utskirts of local galax- 
ies as measured bv lBigiel et al.1 (1201 Obi) . 

On the other hand, the simulations by iGnedin & Kravtsovl 
(2010) focus on galaxies with low metallicities similar to 
those observed of LBGs (~ 0.26 Z Q ) and find reduced SFR ef- 
ficiencies in both their molecular-dominated and their atomic- 
dominated gas. They also find that the efficiency is directly 
related to the metallicity of the gas, suggesting that the de- 



creased star formation efficiency is most likely due to de- 
creased metallicities. These models accurately predict the 
SFR efficiencies we measure in the outskirts of the LBGs 
if they are associated with DLAs. The outskirts of the local 
galaxies measured by Bi giel et ail (1201 Obi) are also generally 
of lower metallicitie s (e.g. Ide Paz et al.1 120071: iCionil 120091 : 
iBresolin et al.1 12009). so the decreased efficiencies could also 
be due to the lowe r metallicity. The SFR efficiencies in 
Bigiel et al. (2010b) are similar to this study, and we discuss 

this below in Section 7.3.1 . 

In addition, iBigiel etafl (I2010al) find a clear correlation of 
the FUV light (representing star formation) with the location 
of the H I gas in M83. In the inner region they find a much 
steeper radial decline in Esfr than in Ee n while in the outer 
region o f the disk, E.sfr decl ines less steeply, if at all (see Fig- 
ure 4 of IBigiel et aLll2010a1) . Similarly, the metallicity in the 
inner part o f M83 drops pretty steeply and then flattens out at 
larger radii (Bresolin et al. 2009). These results are consistent 
with the scenario that the metallicity is driving the efficiency 
of the SFR. 

While it is not yet clear what causes the decrease in SFR 
efficiency in DLAs, the results suggest it is either due to 
a different KS relation for atomic-dominated gas or due to 
the metallicity of the gas rather than the redshift. Given the 
excellent agreemen t of o ur results with the predictions by 
IGnedin & Kravtsovl d2010l) . and the suggestive results of the 
outer regions of M83, we give extra credence to the metallic- 
ities driving the SFR efficiencies. In fact, these two effects 
may be the same, as the outskirts of galaxies generally prob- 
ably have lower metallicities, and therefore lower SFR effi- 
ciencies. While all the gas at high redshift may not have a 
reduced efficiency, care must be taken when using the KS re- 
lation in cosmological models, as the properties of gas vary 
with redshift, thereby affecting the SFR efficiencies. 

7.3.1. Comparison of the z ~ 3 SFR Efficiency with Local z = 
Outer Disks 

The environment for star formation at large galactic radii 
is significantly different than in the inner regions of star- 
forming galaxies, having lower metallicities and dust abun- 
dances, consisting of more H I than H2 gas, and being spread 
out over large volumes. These environmental factors undoubt- 
edly have effects on the conversion of gas into stars, and per- 
haps the SFR efficiency. Similar to our results at high red- 
shift, the SFR efficiency in the loca l (z = 0) outer disks is also 
less efficient than the KS relation (Bigiel et al. 2010b). Fig- 
ure [12] plots our results of H I dominated gas in the outskirts 
of LBGs at z ~ 3 from F igure [TUl in conjunction with the 
measurements of Bi giel et all d20 1 Obi) of the outskirts of local 
spiral galaxies. These measurements combine data from the 
H I Nearby Galaxy Survey and the GALEX Nearby Galaxy 
Survey to measure both the atomic hydrogen gas surface den- 
sity and the FUV emission tracing the SFR in 17 spiral galax- 
ies. The blue diamonds in Figure [T2l represent t he best esti- 
mate fo r the true relation of Esfr and Shi from Bigi el et ail 
(2010b) accounting for their sensitivity, and the error bars 
represen t the scatter i n the m easurements. The SFRs deter- 
mined in Bigiel et al.l d20 1 0b) from the FUV lumi nosity use a 
Kroup a-type initial mass function (IMF) and the Sali m et al.l 
( 120071) FUV-SFR calibration, whil e we use a Salpeter IMF 
and the FUV-SFR calibration from lMadau et all d 19981) and 
iKennic utt (19983). For comparison with ou r LBG results, 
we convert the SFRs from Bi giel et all {2010b) to Sal peter by 
multiplying by a factor of 1.59 and to the iKennicuttl ([l998b) 
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Figure 12. St ar formation rate pe r unit area (Sjpr) versus gas density (S eas ) for the data as shown in Figure [To] and the res ults from the outer r egions of spiral 
galaxies from I Bigiel et al. 2010b). The blue diamonds represent the best estimate for the true relation of Esfr and Shi from Bigiel et al. (2010b) accounting for 
their sensitivity, and the error bars represent the scatter. These points h ave been converte d to use the same IMF and FUV-SFR conversion as used in this work. 
The lo ng d ashed line is the Kennicutt relation for local (z ~ 0) galaxies I Kennicutt 1998a). The gray and gold lines and the black points are the data as shown in 
Figure [lO] The dotted red line represents a power-law fit to the data with both the normalization and the slope as free parameters, and the green dot-dashed line 
is a power-law fit with the slope set to f) = 1.4, as described in Section 7.3.1. 



FUV-SFR calibration by multiplying by a factor of 1.30. 

While these populations probe very different redshifts and 
gas densities, they both probe diffuse atomic-dominated gas in 
the outskirts of galaxies and have low metallicities. In fact, the 
differing redshifts enable a comparison of drastically different 
gas surface densities not possible at low redshift alone, due to 
the saturation of H I gas above a thresh old surface d e nsity as 
described in Section 7.2. Specifically, Big iel et al.l (|2010b) 
cannot probe significantly higher gas densities, as at z = the 
H I gas saturates at threshold surface densities of ~ 10M Q 
pc" 2 (see Figure 8 in Bigiel et al. 2008). On the other hand, as 
described in Section 7.2, at z ~ 3 the H I gas does not saturate 
until much higher gas surface densities. However, our mea- 
surements at z ~ 3 are also limited in the other direction, as 
we cannot probe lower gas surface densities due to insufficient 
surface-brightness sensitivity, and therefore we are left with 
non-overlapping results for 10M Q pc 2 < E gas < 30M Q pc 2 . 

It is not likely a coincidence that the SFR is less efficient in 
the outer regions of galaxies at both low and high redshifts. 
It may be that the SFR is less efficient in the outer regions of 
galaxies at all redshifts, and this possibility needs to be inves- 
tigated. There is evidence that the metallicities in the inner 
regions of galaxies observed in emission are higher than that 
in the oute r parts sampled by quas ar absorption line systems 
at z ~ 0.6 (IChen et alJl2005lE007b . If the SFR efficiency de- 
pends on the metallicity, then we would expect lower efficien- 
cies in the outer parts of these z ~ 0.6 galaxies as observed. 

Comparing the data at low redshift with those at high red- 
shift, we find that they appear to fall on a straight line in log 
space. We therefore fit a power law to the outer disk mea- 
surements in Figure Q~2] using the KS relation in Equation 
(O and setting E c = 1 M Q pc" 2 . Leaving both the normal- 
ization and the slope as free parameters, our least-squares 
solution results in a normalization of K = (1.4 ±0.7) x 10~ 5 
Mq yr" 1 kpc" 2 and a slope of /3=1.3±0.1, and is plotted as a 
dotted red line in Figure [12] In addition, we also fit a power 



law whe re we set the slop e to the same value as in the KS 
relation (lKennicuttlll998allbh . namely ,3=1.4. This fit has a 
normalization of K = (8.5 ± 0.7) x 10" 6 M© yr" 1 kpc" 2 and is 
plotted as a dot-dashed green line in FigureQ~2] The fit is con- 
sistent for both the low- and high-redshift galaxy outskirts, 
and the normalization is significantly lower than that of the 
KS relation, with a K value that is 0.03 times the local value. 
We note that the slope of the LBG outskirts is larger than the 
slope for the outskirts of local galaxies. The metallicities of 
the outskirts of local galaxies are probably somewhat higher 
than those in the outskirts of LBGs, assuming they consist 
of DLA gas, and therefore this difference may be a result of 
differing metallicities. 

We acknowledge that we may be comparing different en- 
vironments when considering the high and low redshift outer 
disks, but they may be similar enough to tell us about star for- 
mation in the outskirts of galaxies. Indeed, both populations 
exhibit similar SFR efficiencies. Measurements of the SFR 
efficiency across a range of redshifts are needed to further ex- 
plore if variations of the slope are real and evolve over time. 

7.4. Bimodality of DLAs 

There is evidence to suggest that there are two populations 
of DLAs (high cool and low cool) , since the d istribution of 
cooling rates of DLAs is bimodal (Wol fe et al.l l2008). These 
two populations have significant differences in their cooling 
rates (and therefore SFRs), velocity widths, metallicity, dust- 
to-gas ratio, and Si II equivalent widths. The gas in the high 
cool population cannot be heated by b ackground a nd/or in 
situ star formation alone, and so I Wolfe etafl (120081) suggest 
that the high cool population is associated with compact star- 
forming bulges, or LBGs. These high cool DLAs typically 
have metallicities (Z = 0.09 ± 0.03 Z^), higher than the av- 
erage for DLAs (IWolfe et al.ll2.Q08b . The low cool popula- 
tion may also be associated with LBGs, although it is possi- 
ble that this population is heated solely by in situ star forma- 
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tion. These low cool DLAs gener ally have lower tha n average 
metallicities (Z = 0.02 ±0.01 Z ? : lWolfe etaT]|2008l) . 

Given that our measurements in the outskirts of LBGs are 
consistent with both metallicities, it is possible that all DLAs 
are in the outskirts of LBGs. Alternatively, it would also be 
logical to associate the high cool population of DLAs with 
the gas responsible for the star formation in the outskirts of 
the LBGs, and the low cool population woul d then naturally 
be rep resented by the DLA upper limits from Wolfe & Chen 
(2006). While those DLAs are also likely associated with 
galaxies, their low er surface brigh t ness m ay result in this be- 
ing included in the Wolfe & Chen (2006) sample. 

If we are mainly probing the high cool DLAs with the ex- 
tended LSB emission around LBGs, and are mainly prob- 
ing the low cool DLAs in the upper limits from Wol fe et all 
(2008), then the number count statistics would need to be 
modified to account for the different percentages of high cool 
and low cool DLAs. However, these percentages are not yet 
very well constrained. A reasonable first try to modify the 
DLA models is to use the current observational result that 
about half the DLAs are high cool and half are low cool. 
If only high cool DLAs are included, it would lower the ex- 
pected //*, which would increase the efficiencies of the SFR, 
moving the data to the left in Figure [TT] Assuming that each 
of the two DLA populations represents half the DLAs, we 
find the data move to the left by about ~ 0.1 dex. This does 
slightly worsen the agreement of the outskirts of LBGs with 
the blue low metallicity model, but not by much. Also, in this 
scenario, the gas would have higher metallicities (Z ~ 0.09), 
and therefore we would expect the outskirts of the LBGs to 
fall to the left of the blue line. Regardless, the uncertainties 
in the models and the data are large enough that they could 
be compatible with a wide range of interpretations. In short, 
our SFR efficiency results do not drastically change whether 
or not there is a bimodality in the DLA population. 

7.5. Effects of Different DLA Sizes 

One caveat remains, which is the possibility that ;/ DLAs 
are not associated with LBGs, and the FWHM linear di- 
ameters of DLAs at z ~ 3 are less than 1.9 kpc, then 
we have not put a limit on the in situ sta r formation in 
DLAs . However, as discussed in §3.1 of IWolfe & Chenl 
(2006), all models suggested so far predict that the bulk 
of DLAs will have diameters larger than 1.9 kpc (e.g., 
I Prochaska & Wolfd 119971: lMoetal.1 119981: lHaehnelt et al l 
2000t iBoissier et all fcOOll N a gamine et all I2006L 120071 
Tescar i et al.ll2009t iHong et al]l2010l), which is in agreemen t 
with the obs e rvatio ns (IWolfe et al.ll2005b ICooke etai]l2010t) . 
iRauch et~aT] (120081) find spatially extended Ly-a emission 
from objects with radii up to 4" (~30kpc) and mean radii of 
~ 1" (~8 kpc) which they argue corresponds to the H I diam- 
eters of DLAs. While continuum emission has not been de- 
tected on these scales, it may have been observed from associ- 
ated compact objects. In fact, almost all of the Ly-a selected 
objects have at least one plausible continuum counterpart (M. 
Rauch, 2010, private communication). 

If we assume that the sources of Ly-a photons are dis- 
tributed thr oughout the gas, th en we can compare the Ssfr 
values from Ra uch et al.1 (120081) to ours. In this scenario, the 
in situ star formation in the DLA gas is the source of the 
Ly-a photons, which is possibly associated with LBGs de- 
tected in the continuum. We note that currently there is no 
evidence for this assump t ion, a nd it is counter to the inter- 
pretation by IRauch et all (120081) . who assume that Ly-a ra- 



diation originates in a compact object embedded in the more 
extended DLA gas. On the other hand, we cannot rule out this 
possibility and it is intriguing, so we consider it here. Using 
the median size of their measured extended Lya emission of 
- 1", we convert their SFRs to S SF r of 4x 10" 4 - 8 x 10" 3 M 
yr" 1 kpc" 2 . This range in Ssfr overlaps our measured Ssfr 
shown in Figures l4l and [TOl suggesting that we may be mea- 
suring the same star formation in DLAs in two very different 
methods. If this is the case, then the sizes of DLAs may be 
even larger than otherwise expected. Whether or not this in- 
terpretation is correct, we are defin itely above the minimum 
size probed by Wolf e & Chenl |2006) of 1 .9 kpc, implying that 
the above mentioned caveat is unimportant. We are therefore 
confident that we have shown in Section 5 and Section 6 that 
the SFR efficiency in diffuse atomic-dominated gas at z ~ 3 is 
less efficient than for local galaxies. 

7.6. The "Missing Metals" Problem 

The metal production by LBGs can be compared to the met- 
als observed in DLAs. Previously, the metal content produced 
in LBGs was found to be significantly larger than that ob- 
served in DLAs by a factor of 10 and was call e d the "Miss- 
ing Metals" problem for DLA s dPettinil 119991 12004 120061: 
lPage| 120021: IWolfe et alJl2003at iBouche et alj(2005l) . How- 
ever, [Wolfe & Chen ( 200q) found that when taking into ac- 
count a reduced efficiency of star formation as found in DLAs 
at z ~ 3, then the metal over production by a factor of 10 
changed to one of underproduction by a factor of three. 

In this study, we calculated the metal production in the out- 
skirts of LBGs and found that the metals produced there yield 
metallicities in the range of 0.12±0.05Z Q to 0.19±0.07Z Q , 
depending on the outer radius for the integration of the SFR 
(see Section 5.3 and Table 13. The average m etallicity of 
DLAs at z ~ 3 is Z ~ O.O4Z (IProchaska et al .1120031 ; M. 
Rafelski et al. 2011, in preparation), and the metallicity of 
the high cool DLAs likely associated with LBGs (see Sec- 
tion 7.4) typically has metallicities of Z = 0.09±0.03 Z . We 
therefore find that if all the metallicity enrichment of DLAs 
were mainly due to in situ star formation in the outskirts of 
LBGs, then there would be no "Missing Metals" problem. 

In order for the metallicity enrichment of DLAs to arise 
primarily from in situ star formation in the outskirts of LBGs, 
we would require some mechanism for the LBG cores, where 
the metallicity is high, to not enrich the intergalactic medium 
(IGM) significantly from z ~ 10 to z ~ 3 (~ 2 billion years). 
However, large scale o utflows of several hundred km s" 1 are 
observed in LBGs (e . g.jFranx et al Jl 1 997h I Steidel et al.ll200l 
Pettini et al.l l2002bt lAdelberger et all 120031: IShaplev et alj 
20031: ISteidel et alj|2010h . Therefore, the lack of enrichment 
would require that either these outflows do not mix signifi- 
cantly with the circ umgalactic medium (CGM, as defined by 
ISteidel et all (120101) to be within 300 kpc of the galaxies) in 
the given time frame or that the outflows do not move suffi- 
cient amounts of metal-enriched gas to the outer regions of 
the LBGs to significantly affect the metallicity. We note that 
we also expect some atomic gas to be in the inner regions 
of LBGs, and therefore a subsample of DLA gas may be en- 
riched locally, although with a smaller covering fraction. 

While there are many theoretical models and numeri- 
cal simulations to understand the nature of DLAs (e.g. , 
l Haehnelt et alJI 1 998t iGardner et al.ll20"0lt iMaller et all 1200 It 
Razoumov et al. 2006; Nagamine et al.l 120071 : iPontzen et alj 
2008; Hon g et al]|2010l) . none of them are able to match the 
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column densities, metallicity range, kinematic properties, and 
DLA cross sections in a c o smo logical simulation. In addition, 
other than Pon tzen et al.l d2008l) . th ey do not reproduce the 
distribution of metallicities in DLAs. iPontzen et alJ (120081) do 
not yet address the mixing of metals nor the effects of large 
outflows on these metals, and can therefore not address this 
question either. These are essential for understanding the pro- 
duction and mixing of the metals, and it is therefore unclear 
at this time if outflows from LBGs mix significantly with the 
CGM, or if they move enough metal content from the inner 
core to the outskirts to significantly increase the metallicity. If 
the outskirts are not significantly enriched by the cores, then 
the observed metallicity of DLAs would be consistent with 
expectations from the in situ star formation in the outskirts of 
LBGs. 

8. SUMMARY AND CONCLUDING REMARKS 

In this work we aim to unify two pictures of the high red- 
shift universe: absorption line systems such as DLAs that pro- 
vide the fuel for star formation, and compact star-forming re- 
gions such as LBGs which form the majority of stars. Each 
population provides valuable but indepen dent information 
about the early universe (for reviews, see Giavaliscol 120021; 
IWolfe et al. 2005). Connecting these two populations helps 
us to better understand how stars form from gas, an important 
part in understanding galaxy formation and evolution. In do- 
ing so, we begin to bridge two separate but complementary 
fi elds in astrophysics. 

Wolf e & Chenl (120061) start to bridge these fields by setting 
sensitive upper limits on star formation in DLAs without 
compact star-forming regions, finding that the in situ star 
formation in DLAs is less than 5% of what is expected 
from the KS relation. However, they do not constrain DLAs 
associated with bright star-forming cores such as LBGs. 
In the present paper we address this caveat by search- 
ing for spatially-extended star formation in the outskirts of 
LBGs at z ~ 3 on scales up to ~ 10 kpc. We find the following. 

1. Using the sample of 407 z ~ 3 LBGs in the UDF from 
iRafelski et al.l (12009ft . we create a composite image stack in 
the V band, corresponding to the rest-frame FUV emission 
which is a sensitive measure of the SFR, for 48 resolved, 
compact, symmetric, and isolated LBGs at z ~ 3 (Figure [3). 
We detect spatially extended LSB emission in the outskirts 
of LBGs, as shown in the radial surface brightness profile 
in Figure [4] This is evidence for the presence of in situ star 
formation in gas in the outskirts of LBGs. 

2. We find that the area covered by DLAs is larger than 
the area of the outskirts of LBGs for SFR efficiencies 
of K = K^enn and is consistent with DLAs having SFR 
efficiencies of K ~ 0.1 x ^Kenn (Figure |5). On the other 
hand, the covering fraction of molecular gas is inadequate to 
explain the star formation in the outskirts of LBGs (Figure 
[6j. This suggests that the outskirts of LBGs consist of 
atomic-dominated gas, supporting the underlying hypothesis 
of this paper. In fact, the covering fraction provides evidence 
that the SFR efficiency of atomic-dominated gas at z ~ 3 is 
on the order of 10 times lower than for local galaxies. 

3. We integrate the rest-frame FUV emission in the outskirts 
of LBGs and find that the average SFR is ~ 0.1 M Q yr -1 and 
p\ is - 3 x 10" 3 M Q yr" 1 Mpc" 3 (see Table |2). Combining 



our largest possible value of p'» in the out s kirts o f LBGs with 
the upper limit found in IWolfe & Chenl (|2006), we obtain 
a conservative upper limit on the total p\ contributed by 
atomic-dominated gas of p'*< 9.9 x 10~ 3 Mq yr" 1 Mpc -3 . 
This corresponds to ~ 10 % of the p % meas ured in the inner 
regions of LBGs at z ~ 3 (|Reddy et alJl2008h . 

4. We integrate p* in the redshift range 3< z <10, and calcu- 
late the total metal production in DLAs and get a metallicity 
of ~ 0.15Z Q . This is comparable to the metallicity of the 
high cool DLAs believed to be associated with LBGs with 
metallicities of ~ 0.09Z Q . If the large observed outflows 
of several hundred km s" 1 of LBGs do not significantly 
contaminate the CGM, then the metallicity of DLAs would 
be consistent with expectations from the in situ star formation 
in the outskirts of LBGs. This is a potential solution to the 
"Missing Metals" problem. 

5. Under the hypothesis that the observed FUV emission 
in the outskirts of LBGs in from in situ star formation in 
atomic-dominated gas, we develop a theoretical framework 
connecting the emission observed around LBGs to the ex- 
pected emission from DLAs. Such a framework is necessary 
to interpret the spatially extended star formation around 
LBGs. This framework develops a differential expression 
for the comoving SFR density, corresponding to a given 
interval of surface brightness using the KS relation and 
the column-density distribution function (Equation (f25ll). 
We also develop a method to convert the measured radial 
surface brightness profile of the LBG com pos ite into the 
same differential expression for p\ (Equation (128b'). 

6. We compare the predictions for the surface brightnesses 
and dp'*/d{I^ s ) to the measured values in Figure [9] and 
find that the two overlap if the efficiency of star formation 
in neutral atomic-dominated gas is lower than the local KS 
relation by factors of 10-50. Using these re duced efficien- 
cies, w e convert our results and those from Wolfe & Chenl 
(2006) into the standard Esfr versus E gas plot generally 
used to study SFRs (Figure ITOb. We find that Esfr is lower 
than the upper limits found by Wol fe & Chenl (120061) for 
a similar range in S gas , and both results have significantly 
lower SFR efficiencies than predicted by the local KS relation. 

7. The reduced SFR efficiencies i n the outskirts of LBGs are 
consistent with the predictions by Gnedin & Kra vtsovl ( 120101) 
for star formation at z ~ 3 in neutral atomic-dominated gas 
with low (< O.1Z ) metallicities (right panel, Figure fill). 
These models find that the primary cause for the lower SFRs 
is due to the decreased metallicities, suggesting that this is 
likely the primary driver for the reduced SFR efficiencies. 

8. Our results correspond to gas surface densities below the 
predicted threshold for saturation of atomic gas at z = 3 for 
DLA metallicities, and therefore the low measured SFR effi- 
ciencies are not due to the saturation of the atomic-dominated 
gas. In fact, our data support the theoretical pr edictions 
dKrumholz etafl l2009at iGnedin & Kravtsovl 120101) that the 
atomic-dominated gas does not saturate for Ehi of at least 
> 1OOM pc" 2 . 

9. Our finding of star formation in atomic-dominated gas 
in the outskirts of LBGs is similar to the recent results by 
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Bigiel et al. (2010b) who find that star formation in the out- 
skirts of local (z = 0) galaxies also arises in atomic-dominated 
gas (Figure[T2"li. In fact, the two results are consistent with the 
same power law, and both results find that the SFR efficiency 
in this gas is lower than expected from the KS relation. It is 
possible that the SFR efficiencies in the outskirts of galaxies 
may be lower at all redshifts, and this tantalizing possibility 
should be investigated. 

10. We find that the reduced efficiencies of star forma- 
tion are likely due to either a different KS relation for 
atomic-dominated gas or due to the metallicity of the 
gas. It is possible that the lower efficiencies are due to 
atomic-dominated gas following a different KS relation than 
molecular-dominated gas, however, we favor the idea that the 
metallic ity of the gas drive s the S FR (e. g., iKrumholz et alJ 
I2009at iGnedin & Kravtsovl l2010l l201lh . In fact, these 
two effects may be the same, as the outskirts of galaxies 
generally probably have lower metallicities, and therefore 
lower SFR efficiencies. While the reduced efficiency may not 
be observed for all gas at high redshift, care must be taken 
when applying the KS relation, as the properties of gas vary 
with redshift, thereby affecting the SFR efficiencies. 

Moreover, simulations of galaxy formation incorporating 
star formation need to take these latest results into account, 
as the local KS relation may not be valid when dealing with 
either atomic-dominated gas or gas with lower metallicities. 
Further observations of star formation in atomic-dominated 
gas are needed to differentiate between these two possibilities 
responsible for reducing the SFR efficiencies. Future stud- 
ies similar to this one at a range of redshifts will help give 
a clearer picture, as will further studies of star formation in 
the outskirts of local galaxies. In addition, progress will also 
hopefully be made by measurements of the C II 158/im line 
emissi on in DLAs using th e Atacama Large Millimeter Array 
(e.g.. Nagamin eetal]|2006h . 
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APPENDIX 

A. COMPLETENESS CORRECTIONS 

Throughout the paper, we use the number of z ~ 3 LBGs in 
the UDF, A^lbg- However, the number of LBGs detected de- 
pends on the depth of our images, and there are large numbers 
of LBGs fainter than our detection limit that are not identified. 
These LBGs are part of the underlying sample and contribute 
to Ca, p*, and Ap*/A(/* s ). We describe our completeness 
corrections for each of these in the following sections. 



A. 1 . Covering Fraction Completeness Correction 

The covering fraction of LBGs depends on Nlbg, and we 
recover the underlying covering fraction for all z ~ 3 LBGs 
by applying a completeness correction which assumes that the 
undetected LBGs have a similar radial surface brightness pro- 
file as the brighter detected LBGs. We acknowledge that the 
size of LBGs vary with brightness, but do not have measure- 
ments of its variation at these faint magnitudes, and there- 
fore make this assumption. If the fainter LBGs were smaller, 
then their \iy would be larger at smaller radii, and it would 
decrease the slope of the covering fraction completeness cor- 
rection. This would not make much of a difference for the 
correction in Figure [5] but it could potentially have a small 
effect on the completeness correction in Figure[6] 

The first step in calculating the completeness correction is 
determining the total number of LBGs expected in each half- 
magnitude bin from the nu mber counts derived fro m the best- 
fit Schechter function from lReddv & Steidell (120091) in Section 
3.1. Since the number counts in Figured] ag ree well with the 
number counts of the sample from Rafelski et al. (200§), we 
are confident that this is the appropriate number of expected 
LBGs. We then subtract the number of detected LBGs in each 
bin yielding the number of missed LBGs per half-magnitude 
bin, Nj(mi S ), where j is the index in the sum in Equation (f30b 
representing the half-magnitude bins. We repeat this for 20 
bins from V ^27 to V ^37 mag. The results are insensitive to 
the faint limit, as described below. 

We then determine the ratio of the flux that needs a cor- 
rection for each half-magnitude bin, F (mag)y, to the total 
flux of the composite LBG, F comp , such that the flux ratio is 
Rj = F '(mag) j/F comp . For each half-magnitude bin, we scale 
the profile by Rj and then determine the area, Aj, that the 
scaled profile covers for each fly. The missed covering frac- 
tion, (CA) m i s , as a function of \iy is just the sum over all half- 
magnitude bins, namely, 

20 

(CUU = J2 a j n j(™-> ■ (30) 

7=1 

We treat the completeness correction differently for atomic- 
dominated and molecular dominated gas. For the atomic- 
dominated gas, we only consider the light from the outskirts 
of these missed galaxies, and not the inner cores. On the 
other hand, for the molecular-dominated gas, we consider 
light from the entire LBG, making no distinction between the 
outskirts and the inner parts of the LBGs. In this scenario, 
both the cores and the outskirts are composed of molecular 
gas, so the cores also contribute to the surface brightness in 
the outskirts for each /iy. The gold long-dashed lines in Fig- 
ures [5] and [6] correspond to the covering fraction for the LBGs 
corrected for completeness. 

We note that while A^(mis) increases steeply with increas- 
ing magnitude, Rjiia&o) decreases more steeply, and therefore 
the completeness correction is dominated by the bright end. 
To give an idea of the magnitude range of the LBGs con- 
tributing the most to the correction, we discuss the correc- 
tion for the entire LBG profile here. The missed LBGs fainter 
than V ~ 33 have basically no effect, since the surface bright- 
nesses barely overlap the area of interest. We could have trun- 
cated the completeness correction at this magnitude with no 
changes to our results. The amount of correction depends on 
the nv being considered, and half the completeness correction 
for the faintest jiy values is due to LBGs brighter than V ~ 3 1. 
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Additionally, those V < 31 LBGs only contribute to /iy > 28. 
In other words, for py < 28, we could have truncated the com- 
pleteness correction at V ~ 31 without any change. 

Hence, the luminosity function used at the fainter magni- 
tudes is not crucial, as those points add negligible amounts to 
the correction. However, if there was a significant deviation in 
the faint end of the luminosity function between 27 < V < 3 1 , 
it would affect our completeness corrections when consider- 
ing both the cores and the outskirts, (see Figure |6}. We note, 
however, that even if we only correct for completeness for 
missed LBGs with V < 29, where a drastic change in the lu- 
minosity function is unlikely, an evolution of f* in /(N^) 
of ^40 would be needed to account for the LSB emission in 
the outskirts of LBGs. Therefore, variations in the faint end 
luminosity function are not responsible for the disagreement 
observed in Figure|6] In the case of the atomic-dominated gas 
where we only consider the outskirts, the completeness cor- 
rection is even more dominated by the brightest missed LBGs, 
since the fainter outskirts quickly do not overlap the p,y of in- 
terest, making the total correction very small (see Figure |5}. 

A.2. p\ Completeness Correction 

Our determination of p\ also depends on A^lbg, and we 
use a similar completeness correction to Appendix A.l here. 
We use the same formalism, using the flux ratio Rj and the 
number of missed LBGs, Nj^) as calculated there. We also 
once again assume that the undetected LBGs have a simi- 
lar radial surface brightness profile as the brighter detected 
LBGs. However, this time in addition to summing over all 
half-magnitude bins j, we also sum over all p,y to get the total 
p\ missed due to completeness, /o'* mis . Specifically, we find 

20 

P*mis = J2J2 P * R J N A^)/ N ^C ■ (31) 
7=1 

We note that similar to Appendix A.l, this result is again not 
sensitive to the correction at the faint end, making the uncer- 
tainty of the luminosity function out there unimportant. 

A. 3. AjO, Completeness Corrections 

The determination of Ap'* depends on Nlbg, which again 
depends on the depth of our images. These missed LBGs 
are part of the underlying sample and contribute to the true 
value of A/?'*, but are not included in Figure [8] We recover 
the underlying true Ap* for all z ~ 3 LBGs by applying a 
completeness correction similar to Appendices A.l and A. 2. 
We use the same formalism developed in A. 1 and once again 
assume that the undetected LBGs have a similar radial sur- 
face brightness profile as the brighter detected LBGs. Just 
like in Appendix A.l, we determine the flux ratio, Rj, and the 
number of missed LBGs, A^(mi s ), where the index j represents 
half-magnitude bins for different missed LBGs. For each bin, 
we calculate the missed Ap* / 'A(7J* 8 ) as a function of py and 
then sum over the half-magnitude bins to get the final correc- 
tion to be applied to Ap'»/A(/^ s ). Specifically, we find 

20 

(^>L = 5 a% r ^™ /n ™ ■ (32) 

We note that similar to Appendix B, this result is not sensitive 
to the correction at the faint end, making the uncertainty of the 



luminosity function out there unimportant. The completeness 
corrected py versus Ap'*/A(/* s ) is shown in Figure|9] 

B. COMPARISON OF THE SUBSET AND FULL LBG 
STACKS 

In Section 3.3, we compared the subset and full LBG sam- 
ples and found them to have similar magnitude, color, and 
redshift distributions. This implies that the two samples of 
LBGs have similar SFRs, stellar populations, and SFHs, and 
justifies using only the subset sample in the analysis to avoid 
contamination. In addition, the KS relation works the same on 
galaxies of different morphologies and environments, so we 
do not expect an effect based on sample selection. Nonethe- 
less, for completeness we investigate how the results would 
differ if we had stacked the full sample of LBGs rather than 
the subsample. A stack of the full sample of LBGs would 
include contamination and different morphological areas as 
described in Section 3.3, and therefore these stacks yield an 
upper limit to the star formation occurring in the outskirts of 
LBGs. 

We repeat the analysis for the full sample stack, and find 
that it yields SFR efficiencies slightly higher than discussed in 
Section 6.2 for the subsample stack. Specifically, the Kenni- 
cutt parameter K needs to be reduced by a factor of ~ 9 below 
the local value for the entire outskirt region. The slope of the 
surface brightness versus Ap*/ A(I°^ S ) is similar to the model 
dp'*/d(I°^ s ), and therefore yields an efficiency that does vary 
with radius. This effect is largely due to contamination from 
nearby galaxies, and if we repeat the analysis for half the sam- 
ple including only isolated galaxies, we find a varying SFR 
efficiency with radius. We conclude that the correct sample 
to stack is the subsample stack, but a stack of the full sample 
does not drastically increase our SFR efficiencies. Therefore, 
the result of lower SFR efficiencies is robust and not due to 
any sample selections. 

We also repeat the covering fraction analysis of Section 
5.2 for the full sample stack and half sample stacks men- 
tioned above. In both cases the SFR efficiencies of atomic- 
dominated gas would need to be increased by a factor of 
about two over the results from the subset sample stack to 
be compatible with the covering fraction of DLAs. Since 
these stacks include contamination and morphologically dif- 
ferent galaxies, molecular-dominated gas from parts of some 
of these galaxies will be contributing throughout. Therefore, 
these results are also upper limits to the covering fraction 
of the atomic-dominated gas. Regardless, even with signif- 
icant contamination, this analysis also indicates that our re- 
sults are robust and not due to sample selections. In addition, 
the molecular-dominated covering fraction based on the full 
LBG stack is still insufficient to account for the emission in 
the outskirts of the LBGs. 
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